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Abstract 

We develop a theory of the domain patterns in systems with competing short-range attractive 
interactions and long range repulsive Coulomb interactions. We take an energetic approach, in 
which patterns are considered as critical points of a mean-field free energy functional. Close to 
the microphase separation transition, this functional takes on a universal form, allowing to treat 
a number of diverse physical situations within a unified framework. We use asymptotic analysis 
to study domain patterns with sharp interfaces. We derived an interfacial representation of the 
pattern's free energy which remains valid in the fluctuating system, with a suitable renormaliza- 
tion of the Coulomb interaction's coupling constant. We also derived integrodifferential equations 
describing the stationary domain patterns of arbitrary shapes and their thermodynamic stability, 
coming from the first and second variation of the interfacial free energy. We showed that the length 
scale of a stable domain pattern must obey a certain scaling law with the strength of the Coulomb 
interaction. We analyzed existence and stability of localized (spots, stripes, annuli) and periodic 
(lamellar, hexagonal) patterns in two dimensions. We showed that these patterns are metastable 
in certain ranges of the parameters and that they can undergo morphological instabilities lead- 
ing to the formation of more complex patterns. We discuss nucleation of the domain patterns 
by thermal fluctuations and pattern formation scenarios for various thermal quenches. We argue 
that self-induced disorder is an intrinsic property of the domain patterns in the systems under 
consideration. 
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I. INTRODUCTION 



Pattern formation is a beautiful example of cooperative behavior in complex systems. It 
is most pronounced in open dissipative systems maintained away from thermal equilibrium 
by external fluxes of energy or matter [|1], ^ ^ ||, ^ At the same time, there exists a 
great number of systems interacting with the outside world only through the contact with a 
heat bath which are also capable of pattern formation and self-organization. Typically, these 
systems are characterized by the presence of coexisting phases, or a phase transition which 
is the driving force for the cooperative behavior. Examples of such classical systems include 
ferroelectric and ferromagnetic films, ferrofluids, Langmuir monolayers, various polymer 
systems, etc. (see, for example, [0, H, |, |10|, [ll], |12|])- Among such quantum systems are type- 
I superconductors in the intermediate state, high-temperature superconductors, degenerate 
ferromagnetic semiconductors, etc. (see, for example |jT3|, |15|, p^). 

In systems not far from thermal equilibrium patterns may form as a result of the com- 
petition of interactions operating on different length scales 0. Typically, a short-range 
attractive interaction in the system would favor macroscopic phase separation. The latter, 
however, is counteracted by a long-range repulsive interaction. This is often accompanied 
by a microphase separation transition, that leads to spontaneous formation of patterns in 
the ideally homogeneous systems upon variations of the control parameters. 

An important class of systems with competing interactions are systems in which the 
long-range interaction is Coulombic. The fundamental nature of the Coulomb interaction 
makes this class of systems extremely diverse. These systems include a variety of poly- 



mer systems, such as block copolymers [QQ, O, 17, ^ M, weakly charged polyelectrolyte 



solutions ||22, ^ cross-linked polymer mixtures p3|; amphiphile solutions phase- 
separating ceramic compounds |^ ; systems undergoing reaction-controlled spinodal decom- 



position photostimulated phase transitions g, |Tj, |2^, etc. Some aspects of systems with 
competing interactions are shared by systems far from thermal equilibrium, such as heated 
electron- hole and gas plasma, semiconductor devices p|, |^ , crystal surfaces undergoing laser- 
induced melting , autocatalytic chemical reactions and surface catalytic reactions [|, p9 



Furthermore, a number of quantum systems, such as degenerate magnetic semiconductors 
and high-temperature superconductors which exhibit electronic phase separation can be con- 
sidered as systems with competing Coulomb interactions O, BO]. In addition, the general 
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problem of Wigner crystallization [|T^, ^ , as well as the thermodynamic and glassy prop- 
erties of spin systems frustrated by Coulomb interaction |3^, ^ |3^, can be 
considered from this point of view. 

Here we develop a theory of patterns with sharp interfaces (domains) in systems with 
short-range attractive interactions and long-range repulsive Coulomb interactions. Our 
starting point is a mean-field free energy functional, which has a nonlocal term associated 
with the Coulomb interaction. Specifically, we are interested in the case of a weak Coulomb 
interaction, when domain patterns with sharp interfaces are realized. We view patterns 
as critical points of the free energy functional. Our main tool in the analysis is singular 
perturbation theory based on the strong separation of length scales in the systems under 
consideration. We use the results of our analysis of the domain patterns to study nucleation 
and formation of complex patterns. We also discuss the effect of thermal fluctuations and 
thermodynamic properties of these systems. 

Our paper is organized as follows. In Sec. II, we introduce the general free energy 
functional and its reduction near a "local" critical point, derive the interfacial representation 
of the free energy and develop a renormalization scheme to account for the effect of thermal 
fluctuations. In Sec. Ill we derive the asymptotic equations for the stationary patterns and 
their stability. In Sec. IV we perform a detailed analysis of localized and periodic patterns 
in two dimensions. In Sec. V we discuss nucleation and growth of complex patterns as a 
result of instabilities of simple patterns, and in Sec. VI we draw conclusions. This paper is 



partially based on the author's Ph. D. Thesis [37 



II. SYSTEMS WITH COMPETING INTERACTIONS OF COULOMB TYPE 
A. Free energy functional 

We start by considering the following general mean-field free energy functional 



a 



+- J d''x'g[<P{x)]G{x - x')g[<P{x')]j . (1) 

Here (f){x) is a scalar order parameter, is a double-well potential, G{x — x') is a positive- 
definite long-range kernel, a is a (positive) coupling constant, (^(0) is a monotone function 



that is equal to zero at some = 0, and d is the dimensionahty of space. Here and henceforth 
we use dimensionless units. 

The functional in Eq. (|1]) may be applicable to a variety of systems. Generally, may 
stand for magnetization, density of the charged polymer in a polyelectrolyte solution, volume 
fraction of a block copolymer in a diblock copolymer melt, density of electrons or holes in 
a charge density wave, structural state of a catalytic surface, concentration of a chemical 



species, etc. [|^, [T^, |T8|, ^ ^7} il, |33, |38|. The kernel G{x — x') we are interested in is the 



Coulomb potential, i.e., it satisfies 

-V^G{x-x') = 6^''\x-x'), (2) 

where S^"^^ (x) is the (i-dimensional Dirac delta-function. The physical nature of the Coulomb 
interaction may also significantly vary from system to system: it may arise as a result of 
the actual electrostatic repulsion due to charges associated with the order parameter, it 
may have an entropic origin, as in block copolymers, or it can come from the diffusion of 



chemically reacting species (see, for example, 0, |T^, ^ ^). Note 

that in quantum systems Eq. (|1|) arises within the framework of density functional theory 
|g,^,|D|. 



The long-ranged nature of G{x — x') from Eq. (|^) is expressed in the fact that its Fourier 
transform has a singularity at wave vector k = 0. At the same time, this Fourier transform 
is bounded at large fc-vectors. Let us emphasize that G{x — x') represents a repulsive long- 
range interaction, since it is positive for all x. Therefore, the Coulomb long-range interaction 
represented by G{x — x') is competing with the short-range interactions represented by the 
first two terms in Eq. (|l]). It is also clear that since the Fourier transform of G{x — x') is 
positive for all wave vectors, the functional in Eq. (|l]) is bounded from below on any finite 
domain. 

If we formally put a = in Eq. (|l]), we will recover the standard free energy functional 
that is used in the studies of phase separation (see, for example, ^1|]). On the other hand, 
no matter how small the value of a is, because of the singularity of the Fourier-transform 
of G{x — x') at k = the effect of the long-range interaction will remain significant on 
sufficiently large length scales. Indeed, if the system has a finite size L, from the dimensional 
considerations the contribution of the Coulomb interaction into the free energy will scale as 
aL'^'^'^. If the value of a is decreased while L remains fixed, the contribution of the Coulomb 
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interaction goes away. This means that when a <^ 1, the system behaves locally as if it 
did not have the long-range interaction. On the other hand, for an infinite system this 
interaction is always relevant since its contribution scales as L'^^'^ > L'^. Therefore, for 

a < 1 (3) 

the long-range interaction will be a singular perturbation, globally affecting the behavior of 
the system. It is in this case that domain patterns form in systems with the free energy of 
the form of Eq. (|l]). Since we are interested in the domain patterns here, Eq. (|^) will be 
assumed from now on. Note that this condition is satisfied in many systems with long-range 
interactions of Coulomb type |jl5|, |18|, ^ |2^, |30l . 



The singularity of G{x — x') on the large length scales implies that the Fourier component 
of g{(f)) at k = must vanish in order for the last integral in Eq. (P to remain finite. This 
corresponds to the overall electroneutrality for systems in which the order parameter is 
associated with the electric charge. The only possible homogeneous phase of the system is, 
therefore, = 0. Thus, due to the long-range interaction the global phase separation in the 
system becomes impossible. On the other hand, as we will see below, the system described 
by the free energy functional from Eq. (0) may be in a patterned state. By patterned states 
(more precisely, by stationary patterns), we will mean the inhomogeneous distributions of 
the order parameter which are critical points of the functional F. 

B. The microphase separation transition 

Let us assume that in the absence of the long-range interaction the system would possess 
a critical point at temperature T = Tc. Then, in the vicinity of Tc the function /(0) can be 
expanded as 

— + (4) 



where t = {T — T^jT^ is the reduced temperature, and a and h are positive constants Wl 



In the following, we will talk about as the "local" critical temperature. Near T^, the value 
of 101 ~ 00 = ('^kl/^)^^^ ^ 1 If also 101 <^ 1, we can expand the function g[<\)) in 



a Taylor series and retain only the first term, so (7(0) ~ const x (0 — 0). Then, rescaling 
the order parameter and length with the values of 0o and the short-range correlation length 
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^ = |ar| 1^^, we can write the free energy from Eq. (|I]) below in the following 



universal form: 



V 2 2^4 



d 



+- J dV(0(x) - - x'mx') - <p)j , (5) 

where we absorbed a constant factor into the definition of F. Here the parameter e, which 
plays the role of the effective coupling constant of the long-range interaction, is given by 

m 



ar 



a'/'\T\-\ (6) 



Notice that in Eq. the value of has been rescaled as well, so it now depends on 
temperature: 

Also, as was discussed above, for Eq. the singularity of G{x — x') at small wave vectors 
implies that the total amount of the order parameter must be conserved (the "electroneu- 
trality" condition): 



1 

V 



>d'^x = (f), (8) 



where V is the system's volume. 

Let us consider small fluctuations of the order parameter 5(f) = cf) — cj) away from the 
homogeneous phase for T < T^.. From the second variation of F from Eq. (|^), the Fourier 
transform of the pair correlation function of such fluctuations is 

<l**^l^> |kp + 3^^-l+e^|k|-^ - t**) 
This correlation function has a maximum at non-zero /c- vectors with |k| = kc-, where 

kc = (10) 

The fluctuations at k^ diverge when = ±|0c|, where 
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The divergence of the fluctuations at |k| = fee signifies an instabihty of the homogeneous 
phase and leads to the microphase separation |^. Note that the instability can only be 
realized if e is small enough; in terms of temperature, it occurs at T slightly below when 
Eq. (D holds, see Eq. (§. 

As the temperature is decreased, the value of e gets smaller. Note that for small a one 
can still be close to Tc and yet have e <^ 1. In this situation the long-range interaction 
can be a singular perturbation (in the sense discussed earlier) even in the vicinity of the 
transition. As was already mentioned, this is a necessary condition for the existence of the 
domain patterns, so below we will concentrate on the case e -C 1. For e <^ 1 the instability 
of the homogeneous phase occurs close to the classical spinodal of the Ginzburg-Landau free 
energy: ~ 1/V3, see Eq. (j^. In this case, according to Eq. (p!0D , the instability occurs 
at kc -C 1. 

There are two regions in the fc-space in which the fluctuations of the order parameter 
around the homogeneous phase, with |0| > |0c|, behave differently when e ^ 1. According 
to Eq. (^, for |k| ~ 1 one could neglect the long-range contribution, so the fluctuations 
(|'50kP) oc V^/(|kp + m^), where = 30^ — 1, are those of the (mean-field) critical phe- 
nomena with the length scale independent of e: 

/ ~ 1. (12) 

On the other hand, for |k| ^ 1 one can neglect the |kp term, so the fluctuations behave like 
{\S(j)\l} oc F/(m^ + e^|k|~^) = iV/m^) [1 — e^/(e^ + m^|kp)]. The first term in this expression 
represents local order parameter fluctuations, while the second is the familiar Debye-Hiickel 



correlation function The length scale associated with the latter is the screening length 

L ~ e-\ (13) 

For e -C 1 the (generally, metastable) equilibrium state of the system should be a domain 
pattern made up of domains of large size ~ R separated by narrow domain walls of width 
~ /. Clearly, the long-range interaction cannot significantly affect the local profiles of the 
order parameter; however, it can affect the locations of the domain walls. The size of the 
domains will be determined by the competition between the surface energy of the domain 
walls ~ i?'^"^ per droplet and the energy of the long-range interaction ~ e^R^'^^, so the 
characteristic size of the equilibrium domain pattern will be of order (in the context of block 



copolymers, see also |1T8| ) 

R ~ e-2/3. (14) 
Note that this result for the global minimizers of the sharp interface limit of Eq. (^) was 



recently proved by Choksi |53]. Choksi also obtained rigorous upper and lower bounds on 
the energy of global minimizers of Eq. in the situation when the screening effects are 
neghgible. 



According to Eq. (pUj), the wavelength of the fluctuations with respect to which the 
instability of the homogeneous phase is realized is 

A = 27T/kc ~ e"i/l (15) 

Comparing all these length scales, one can see that for e <^ 1 the following hierarchy holds: 

/ < A < i? < L. (16) 

This is a crucial property of systems with weak long-range Coulomb interaction. 

C. Interfacial representation of the free energy 

The solutions of the Euler-Lagrange equation obtained from Eq. may be analyzed by 
singular perturbations theory in the asymptotic limit e 0. We will perform this analysis 
in Sec. |111 A| . Now, however, we will use a different method which gives the free energy of 
the domain pattern in terms of the locations of the domain interfaces [^]. This method 
was used by Goldstein, Muraki, and Petrich for a reaction-diffusion system with a weak 



activator- inhibitor coupling Here we develop a procedure that allows to calculate 

the free energy of a domain pattern which takes into account the screening effects. 

Because of the strong separation of length scales we can introduce the following ansatz 
for the distribution of the order parameter: 

(Pix) = (pshix) + (psmix), (17) 

where 0sh represents the sharp distributions, whose characteristic length of variation is of the 
order of the domain wall width (which in our units is of order one), and 0sm represents the 
smooth distributions, whose characteristic length of variation is comparable to the domain 
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size R. The distribution 0sh is chosen in such a way that it is equal to +1 inside the positive 
domains and —1 outside, whereas at the interfaces it is close to the one- dimensional domain 
wall of the Ginzburg-Landau theory [ ^ : 

0sh = tanh-^, (18) 
v2 

where p is the distance from a given point to the interface, which is positive (negative) in 
the positive (negative) domains, respectively. Thus, the location of the interface is built 
into the definition of 0sh- The contribution from to the free energy, coming from the 
integration in Eq. (^ in the vicinity (of order 1) of the interfaces, gives the surface energy 

dS, (To = . (19) 

Here the surface integral gives the total surface area of the domain interfaces and (Tq is the 
surface tension coefficient of the domain wall in the Ginzburg-Landau theory . 



To find the smooth distributions 0sm away from the interfaces, we minimize the free 
energy in these regions. Taking into account that 0sm varies slowly on the length scale of 
order 1, we can neglect the term arising in the Euler-Lagrange equation and obtain 

+ - 03 - y d'^x'G{x - x') (0(x') - 0) = 0, (20) 

where /i is the chemical potential (a constant) coming from the constraint given by Eq. (^. 
On the scale of the domains, 0sh = ±1 away from the interfaces. We will assume that inside 
the domains |0sm| ^ 1, which is justified for R ^ e^^ (see below). This allows us to linearize 
Eq. ( pOD around (p^h away from the interfaces. Using Eq. (|l^ with (pl-^ = 1, Eq. (|20|) is 
written as 

where we introduced an effective field 

^ = -p + e^ f d'^x'Gix - x')(0sh(a:') + (psU^') - 0)- 



(22) 

Note that the constant is basically the coefficient of linear response for the local theory. 
Applying to Eq. (|2^ ) and using Eqs. (0) and (|21|), we obtain 



V2V; + e2«;2V' = e2(0sh-0). (23) 



Note that our definition of ip, togetfier witli Eq. (pTD, implies tliat Eq. (H) is automatically 
satisfied to the leading order in e. This can be seen by integrating Eq. ( |23D over the volume 
of the system and taking into account that for no-fiux or periodic boundary conditions the 
surface integral in the obtained expression vanishes. Also note that for the same reason /i 
drops out from this equation. 
The solution of Eq. (^) is 

^ = y - x')(0sh(a;') - 0), (24) 

where is the screened Coulomb interaction that satisfies 

-V^G, + e^K'G, = 5^''\x-x'). (25) 

These are explicitly given as follows 

2^ exp(— e^la; — x'l) in d=l, 
G,{x-x')= I l-Ko{eK\x-x'\) in d = 2, (26) 



exp{—eK\x—x'\) 



in d = 3, 



4tt\x—x'\ 

where Ko{x) is the modified Bessel function. Thus, the fluctuations of the order parameter 
in the bulk indeed screen the interaction on the length scale L ~ e^^. This means that the 
finite size effects will become unimportant if the system size is much greater than L. Notice 
that the value of ip is estimated as ip ^ e^R^ <^ 1 for i? ^ e^^, justifying the linearization 
used in the derivation. Also, according to Eqs. ( pTD and (0), in this situation the deviation 
of from ±1 is small away from the interfaces. 

Let us now calculate the contribution from the long-range interaction to the free energy. 
Once again, neglecting the |V0p term, expanding the nonlinearity in Eq. around 0sh up 
to the second order in (psm, and taking into account that (p"^^ = 1 away from the interfaces, 
to the leading order in e we can write the contribution of the long-range interaction (up to 
an overall constant) as follows: 



long— range 



e2 



dVx'iMx) - 0)G,(x - x')(0sh(a:') - 0), 

(27) 
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where we used Eqs. (H), (pT]), (p2D, and (p4D. One can see from this equation that the 
screening represented by 0sm enters the free energy only via Ge{x — x'). The integral in 
Eq. can be transformed to an integral over the domain interfaces by using Eq. ( p5[ ) 
and the fact that (p^h = ±1 in the positive (negative) domains After calculating the 
respective integrals and collecting all the terms in the free energy (see Appendix we 
obtain 

+2e^ [ [ (fx(fx'G,{x - x') (28) 

= aQ ^ dS ^ y dS{x ■ h) 

-^jdsj dS'{n ■ n')G,{x - x'), (29) 
where VL^ denotes the positive domains, n is the outward normal to the interface of 



and the surface integrals are over the interface. The first integral in Eq. (^9]) is the overall 
surface area of the interfaces, the second gives the total volume of the positive domains, 
and the third is the non-local contribution of the screened long-range interaction (note the 
distinction with |^). Thus, Eq. (|29|) gives the free energy of the domain pattern in terms 
of the locations of the interfaces only. Note that the unscreened version of Eq. ([28| ) was 
recently derived rigorously by Ren and Wei in the context of F-convergence [i7[ . 



D. Renormalization 

The treatment above is based on the mean-field free energy functional from Eq. (^) and 
therefore neglects the effects of thermal fluctuations. In a fluctuating theory, this functional 
will become an effective Hamiltonian, in general requiring an appropriate field-theoretic 
treatment. Nevertheless, we propose that the effect of the critical phenomena fluctuations 
can be taken into account by an appropriate renormalization of the main parameters of the 
free energy in the interfacial regime |^ . 



Indeed, if one looks at the singularly perturbed [Eq. (|^)] fluctuating system near Tc, 
on small length scales one will see critical phenomena fluctuations of a second-order phase 
transition without the long-range interaction. This will happen as long as the characteristic 
screening length L of the long-range interaction is much greater than the correlation length 
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^ of the critical phenomena fluctuations. The critical exponents associated with the local 



critical phenomena fluctuations must be those of the dimensional Ising model |^. So, 
the local average of the order parameter will be close to a constant = ±00^1^) where 
r = (T — Tc)/Tc is the reduced temperature and (3 is the respective critical exponent. Also, 
the surface tension coefficient of an interface in which the order parameter changes sign is 
(T = aQ\T\'^^^~^\ where u is another critical exponent, and its width is roughly the correlation 
length e = ^ok|-'^ H. 



Observe that the long-range coupling involves integration over regions of size ~ R which 
for the domain patterns must be much greater than the correlation length. Therefore, it is 
the average value of the order parameter that gives the main contribution to the long-range 
interaction energy for R ^ ^. This energy has to be compared with the surface energy, so 
in equilibrium we obtain 

Rescaling the order parameter and length appropriately and introducing the renormalized 
coupling constant 



e 



2 = a02^3|^|2/3-(<^+2)^ (31) 



we can still write down the interfacial free energy of the system in the form of Eq. (^), where, 
as usual, we dropped the primes and neglected an overall constant factor. Caution, however, 
is necessary here in considering the screening effects. As was noted earlier, the constant k 
appearing in the mean-field definition of the screened long-range interaction G^{x — x') is 
related to the coefficient of linear response for the local theory. When the critical phenomena 
fluctuations are taken into account, the value of k can be calculated via the linear response 
function x = Xo\t\~'^ that relates the unsealed values of (psm and ip (below Tc), see Eqs. ( ^1]) 
and (|2^) . After an appropriate rescaling and using the definition of e from Eq. (|3l|) , we 



obtain that 



which is a constant of order one, independent of temperature. In writing Eq. ( p2l) we used 
the scaling relation 2(3 + '-f — du = between the critical exponents |^ . 



Thus, we can renormalize the parameter e and redefine the parameter n to obtain once 
again the free energy of a domain pattern in the form of Eq. (^9j) even in the case of the 
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system locally experiencing critical phenomena fluctuations. Let us point out that the free 
energy of a domain pattern of size of order R ^ ^ will be much greater than ksTc (in the 
unsealed units), see Eq. (pQ]), since a^'^'^/ksTc ~ 1 [0. This leads us to the conclusion 



that in a strongly fluctuating system the domain patterns should be essentially described 
by the interfacial mean-fleld theory, and all the properties of the domains in the fluctuating 
system will be equivalent to those of the domains in the mean-fleld systems described by 
Eq. @, provided that one uses the renormalization given by Eqs. (^Tf ) and (^). Thus, the 
universality discussed earlier for the mean-field model should in fact extend to all systems 
near the local critical temperature as long as the coupling constant a of the long-range 
interaction is small enough. Note that in this situation Eq. @ may be used as a phase-field 
model representation for the free energy of the domain patterns [Q. 

In the renormalization of the main parameters of the system we made an assumption 
that the size of the domains must be much greater than the correlation length ^. According 
to Eq. (Q), this condition is satisfied as long as e ^ 1. In view of Eq. (^l]), this is the case 



when the reduced temperature r is much lower than r = —Tc, where 

^^_^l/[.(rf+2)-2/3]^ (33) 

at which e ~ 1. When the temperature is decreased below —Tc, the long-range interac- 
tion becomes progressively more and more relevant at long distances, until for r ^ —Tc 
(what means e ^ 1) long-lived domain structures with the free energy cost of each domain 
AF/ksTc ^ 1 will start to form. On the other hand, for r ^ +Tc the long-range coupling, 
which scales as a^^, will be much smaller than the effective local coupling, which is of order 
~ I'^r, so one will observe only the critical phenomena fluctuations above Tc. It is a 
question whether there is a microphase separation transition from the homogeneous to the 
patterned phase (which is analogous to the freezing transition in liquids) or there is a smooth 
crossover from one to another in a strongly fluctuating system. It is clear, however, that the 
uniform phase must be thermodynamically unstable when r < — Tc. At the same time, at 
r ~ — Tc the fluctuations are strong, so one can envisage the system as a collection of domains 
that randomly appear and disappear in different locations and move about as particles in a 
dense liquid. In any case, there must exist a narrow transition region —t^ ^ t ^ Tc, upon 
going through which the phase should change from uniform to patterned. 



13 



III. PROPERTIES OF THE DOMAIN PATTERNS 



A. Equations for stationary patterns 

The stationary patterns in the mean-field model introduced in Sec. [11 B| must satisfy the 



Euler-Lagrange equation obtained from Eq. 



.,-.= /A'G,.-.')(^M-^)^0. (34) 

where is the chemical potential. Formally, this integrodifferential equation can be rewritten 
as a pair of stationary reaction-diffusion equations of activator-inhibitor type 0, ^ ^ . 
Indeed, if the last term in Eq. ( ^i)) is denoted by this equation can be rewritten as 

VV + - 0^ - ^ = 0, (35) 
VV + e'(0-0) = 0, (^) = -/x, (36) 

where (■) denote averaging over the system's volume. 

Reaction-diffusion equations of the type of Eqs. ( pSj) and (p6D have been studied by many 
authors (see, for example, 0, |5y, 0, 0, 0]). In the limit e — > their solutions can be 
treated by the methods of singular perturbation theory (matched asymptotics) 0, 0, ^ 



According to singular perturbation theory, the solution (p of Eq. (^) can be broken up into 
the inner and outer solutions. The inner solution varies on the length scale of order 1 and 
describes the variation of the order parameter in the vicinity of the domain interfaces, while 
the outer solution varies on the length scale R of the order of the characteristic size of 
the domains and describes the variation of the order parameter away from the interfaces. 
Similarly, the solution ip of Eq. (^) will vary on the length scale R. 

Since the variable ip varies slowly on the inner scale, it can be considered as constant in 
the interface. Since the curvature of the domain wall is also much smaller than the domain 
wall width, to the leading order we can write Eq. ( ^Sf ) in the vicinity of the interface as 

p-2H| + ,-#-*^0. (37) 

where p is the distance from a given point to the interface, which is positive if the point is 
inside the positive domain and negative otherwise, H = |(/ci + ^2) is the mean curvature 
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of the interface (positive if the positive domain is convex, ki and ^2 are the principal cur- 
vatures), and ipi is the value of ip on the interface. In the following, we will write all the 
formulas in the three-dimensional case, in two or one dimensions one has to set one or two 
principal curvatures of the interface, respectively, to zero. 

Equation (^) can be solved exactly, its solution has the form 0(p) = atanh6p + c, where 
a, b, and c are certain constants. This solution exists only when 2H — ^-^^ ~ ~'^'ipi/V^- 
Since in the domain pattern if ^ 1, it is sufficient to linearize this equation with respect to 
H, so we obtain (see also [06|) 



aoH = -ijj,, (38) 



where ctq is given by Eq. (19). This, in turn, implies that <^ 1 in order for a pattern to 
be stationary. Note that to the leading order 0(p) in the interface is given by Eq. (plsl). 
Away from the interfaces (on the outer scale) (p varies slowly, so one can neglect the 



gradient square term in Eq. (|35D. Then, according to Eq. (^), we have (f) — (f)^ = ip . Since 
we must have 1, this equation can be linearized with respect to (p around (p = ±1 in 
the positive and negative domains, respectively. So, one obtains that to the leading order 



±1 — n'^ip in the outer regions. Here, as in Sec. [11 (J| , we have = |. If we substitute 



this expression into Eq. (|36|), we will obtain precisely the same equation for ^p as Eq. (p3[), 
with (psh = ±1. 

The solution of Eq. (^) can be written as an integral over the domain interfaces (Ap- 
pendix |^, see also ^): 

^ = + ^<f dS'{n' ■ V'(a - G)}, (39) 

where V is the gradient with respect to x'. Combining this equation with Eq. (|38|), we 
obtain the following equation for the locations of the interfaces 

aoH = ^ - ^ ldS'{n' ■ V'(G, - G)}. (40) 

K J 

This equation can be further simplified if the distance between the points on the interface is 
much smaller than e~^. In this case one can expand in Eq. (PO]) in en\x — x'\ and retain 
the terms up to the second order. It is easy to see that only the terms of the second and 
higher orders of the expansion of G^ in e will give non-trivial contributions to the right-hand 



15 



side of Eq. (^OD. Also, in view of the approximations used to derive Eq. (^OD, this equation 
is vahd when the characteristic size R of the domains satisfies 1 <^ i? ^ e^^. 

Equation ( ^OD describes the pressure balance across the interface. Indeed, the term in 
the left-hand side of this equation is the Laplace law, the first term in the right-hand side 
gives the bulk pressure, and the second one gives the nonlocal contribution to pressure due 
to the interaction of the domain walls with each other. 

Let us emphasize that Eq. (^Of ) can also be straightforwardly obtained by computing the 
first variation of the interfacial free energy given by Eq. (|29| ) (see Appendix Therefore, 



this equation also remains valid in the fiuctuating system considered in Sec. p D| . Also, 
note that since the solutions of Eq. ( ^41) in the form of stationary domain patterns can be 
written in the form of Eq. for e ^ 1, Eq. (|29D gives the asymptotic expression for the 
free energy of these patterns. 

B. Deformations of the domain interfaces 

The solutions of the Euler-Lagrange equation given by Eq. ( |5^ ) are critical points of the 
free energy functional from Eq. (|]). Similarly, the solutions of Eq. (|40|) are critical points of 
the interfacial free energy from Eq. (p9D and correspond to the solutions of Eq. ( ^41) in the 
limit e — s> 0. Both these solutions define (generally, metastable) stationary patterns. The 
question, however, arises as to when these patterns are thermo dynamically stable. Since, 
apart from the nucleation phenomena discussed in Sec. |V A| , the effect of thermal fluctu- 
ations is small in both cases, the thermodynamic stability of the patterns is determined 
by the second variation of the free energy functional. Thus, the thermodynamically stable 
stationary patterns will be local minimizers of the free energy. 

The problem of flnding the second variation of the functional in Eq. (^) reduces to the 
calculation of the spectrum of linearization of Eq. (p^. It is not difficult to see that it 
is equivalent to the problem of linear stability of stationary patterns in systems obeying 
gradient descent dynamics. Such a stability analysis in the context of general reaction- 
diffusion systems of activator-inhibitor type was performed in |5^. Here, instead of analyzing 



the second variation of F from Eq. (^, we will use the interfacial free energy from Eq. ( ^91) 
for flnding the spectrum of the fluctuations of the pattern's interfaces. These are, in turn, 
the lowest-lying modes of the spectrum and, therefore, cost the least free energy in Eq. (^. 
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Both these approaches give the same results in the Umit of small e. 

Let us now proceed with the calculation of the second variation of the interfacial free 
energy from Eq. (|29|). A small perturbation of the domain shape means a slight shift of the 
interface in the normal direction by p{x), where x denotes a point on the interface. In terms 
of p{x), the second variation of the free energy from Eq. is (see Appendix 0) 



+^ j dS p^{fi ■ V) j dS'{n' ■ V'(Ge - G)} 

+4.= /../.S'G.(.-.0.(.)p(.'). (41) 

where V is the gradient in x', K = kik2 is the Gaussian curvature at a given point on the 
unperturbed interface, V_l is the gradient along the interface, and the integration is over 
the unperturbed interfaces. 

Different terms in the integrand of Eq. (|4l| ) represent competing tendencies that stabilize 
or destabilize the patterns. The cro|V_Lpp term coming from the surface tension penalizes the 
distortions of the interfaces; the term involving the curvatures 2ao{K — 2H^)p'^ = — o"o(^i + 
^2)P^ < 0, is a destabilizing term coming from the curvature of the interface; the term from 
the second line in Eq. (^) can be rewritten as 2(n ■ Vip)p'^, where ip is given by Eq. (|39[) 
(see Appendix ^), and represents the change in the free energy due to the motion of the 
interface in the fixed effective field ip, this term should be destabilizing also since we would 
generally expect the gradient of to be directed inward at the interface; and the last term 
is a stabilizing action of the long-range interaction. 

To gauge the relative strengths of these terms and determine whether a pattern is stable, 
we need to solve the following eigenvalue problem obtained from Eq. (PT|): 



Lp = Xp, where 

Lp = -ao^lp + ^aoKp -AaoH'^p (42) 
+2{n ■ Vip)p + 4e2 / dS'G,{x - x')p{x'). 



The spectrum of the operator L for a stable pattern should not contain any negative 
eigenvalues. We will analyze the spectrum of L for simple geometries below (see also 



4g, ^ |53l)- Now, however, let us discuss some general properties of 5 F in 



Eq. (^11). It is easy to see from Eq. (PD that a typical size i? in a stable stationary pattern 
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must have the same scahng as that in Eq. ([T^) (this point was first argued in ^ |57 



based on the stabihty analysis of the localized and periodic patterns). Indeed, suppose a 
pattern is made of a collection of droplets of size and distance between each other of order R 
(here for definiteness we will consider three-dimensional patterns). Let us first assume that 
the droplets are too small, so R <^ e~^/^. Consider a fluctuation that increases uniformly 
the volume of one droplet while decreasing the volume of another next to it, so that the 
net volume change is zero (repumping, see [@, ^). Then, if i? ^ e~^/^, the stabilizing 
contribution ~ e^R^ from the last term in Eq. (41) is negligible compared to the destabi- 
lizing contribution from the curvature terms ~ 1, while the |V^pp term is identically zero. 
Therefore, such a fluctuation will lead to the free energy decrease. 

Now, suppose that the droplets are too big, so R^ e^'^l'^ . Let us now perturb the interface 
of one droplet in a localized fashion in the region of size £ -C i?, once again, maintaining 
the overall volume the same (distortion, p, Q). Then the last term in Eq. (^T[), which 
is ~ e^£^ will once again be negligibly small compared to the term from the second line of 
Eq. (0), which is ~ e^PR. On the other hand, the gradient square term in Eq. (|4T|), which 
is of order ~ 1 will not be able to compensate that contribution, if e'^R^^I'^ <^ £ <^ R. 
Such i can always be found when R ^ e^^^^, so this kind of a fluctuation will lower the free 
energy, too. Note that for R ~ this instability result was also obtained by Nishiura and 
Suzuki 



The arguments above lead to an important conclusion that (perhaps, apart from some 
logarithmic factors, see below) the stable stationary patterns must obey the equilibrium 
scaling from Eq. (0), which was obtained on global energetic grounds. In other words, not 



only the global minimizers of the free energy but all local minimizers must generally 
obey this scaling. Note, however, that these arguments do not apply in one dimension (see 
also p|, |50|, |51|, 0). Similarly, the equihbrium scaling from Eq. ([T^) is not necessarily obeyed 



by all stationary patterns, see for example. Sec. [IV A| and P^V B| , contrary to the statement 
of [||. 



IV. EQUILIBRIUM PATTERNS AND MORPHOLOGICAL INSTABILITIES 

Let us now use the tools developed in the preceding sections to analyze the stationary do- 
main patterns with simple geometries, such as localized and periodic patterns. In this paper, 
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we will limit ourselves to studying two dimensional patterns. Qualitatively the same results 
are expected for the more experimentally relevant three-dimensional patterns. Note, how- 
ever, that the one-dimensional case is qualitatively different from two and three dimensions 
(see 0,0, 

Since our system possesses a symmetry (j) —(p, we can only consider the properties of 
the positive domains immersed in the negative background. This means that we only need 
to study the region of the system's parameters in which < 0. 

A. Solitary patterns 

We begin with the study of the simplest possible domain patterns: solitary patterns. In 
two dimensions we will consider spots, stripes, and the annuli. 

1. Spot 

Let us first look at a spots: a small positive circular domain in the negative background. 
If the radius of the spot is much smaller than the screening length e~^, the interaction 
potential G^{x — x') in Eq. ( p9D can be expanded in e. Retaining the terms up to e^, 
after a straightforward calculation we obtain that the free energy of a spot of radius TZ is 
asymptotically 

Fin) = 27^ao7^ ^ 



-7^e27^Mln(^^e«7^^ +7-^1 , (43) 

where 7 0.5772 is the Euler constant and 

6 = 1 + ^ (44) 

measures the degree of metastability of the homogeneous phase. The free energy of the spot 
given by Eq. (^31) for a particular set of parameters is shown in Fig. |l]. 

Let us now analyze Eq (P3|). First of all, when 6 < 0, the free energy is a monotonically 
increasing function of TZ. When the value of 6 is increased, at 5 = 5^ ^ 1 a minimum and 
a maximum of the free energy appear (see Fig. |1]). These correspond to the radially stable 
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and unstable spot solutions, with the radii IZ = lZs and TZ = Tin, respectively. At S = Sb we 
have TZs = Tin = 'T^min- Asymptotically 

= (^^eK^ao In^/' e''^ ^'^ ~ e^/^ c\ (45) 

3ao ^'1' 



This formula agrees up to the logarithmic factor with Eq. (|T^. These logarithmic factors 
are a specific feature of the two-dimensional patterns, they are absent in three dimensions 
5|. 



When the value of b is increased beyond 5b, the radius Tls grows, while the radius Tin 
shrinks. At some value of 5 = 5^ •'^ 1 at which Tl^ = Tl^, the free energy of the spot becomes 
negative, making the spot thermodynamically more favorable than the homogeneous phase. 
Once again, asymptotically, 

S,„ = i (se^Vo W'^- ~ e-i, (47) 



Comparing Eqs. (HSl) - P3), we see that 6m = 2^/^5b, and 7^„ = 2^/^Tl^i^. 
For 5 ^ 5fe the radii Tig and T^n become asymptotically 

This means that for S ^ e^^^lne^^, the radius Tig goes beyond the equilibrium scaling of 
Eq. (|Tip. This is an indication of a morphological instability studied in Sec. |1V (J| . 

2. Annulus 

Let us now analyze the pattern in the form of a thin annulus, which has the radius Tl and 
thickness C <^Tl. Calculating the free energy of such a pattern from Eq. (^), we obtain 

F(7^, C) = 47^(To7^ - ^TIC 

+ A-nt^Tl^C'^I^iitKTViK^itKTl), (50) 

where Iq{x) is the modified Bessel function. Minimizing this expression with respect to £, 
we obtain that the value of £ = £a in equilibrium is related to Tl as follows 

" 2e^Kmio{eKTl)Ko{eKTl)' ^^^^ 
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Substituting this expression into Eq. (0), we then study the critical points of F with respect 
to 7^. 

The analysis of Eqs. (^) and ( |5TD shows that for 5 <^ e^^^ there exists a single mini- 
mum of the free energy, corresponding to an annulus solution, whose radius and width are 
asymptotically 

7^a = , Ca = ln(er2). (52) 

One can see that the condition <^ TZa used in the derivation of Eq. (|52D is satisfied as 
long as 5 ^ e^^^ In e~^. 

According to Eq. (|5^) , when 5 ~ e^/^, we have TZa ~ e^^, so screening effects become 
important. The analysis of Eq. ( ^OD shows that at some critical value of 5 ~ e^/^, a new 
minimum and a maximum of F{TZ) appear [5 = 0.0255 in Fig. H). At a slightly higher 
value of (5 ~ e^/^, the second minimum of the free energy disappears {6 = 0.0258 in Fig. 
The value of 6 = 6± at which this happens can be easily calculated, see Eq. (0). 



3. Stripe 



Let us now determine the equilibrium parameters of a quasi one-dimensional domain 
pattern — stripe. A stripe of width Cg can be considered as a limit of an annulus as 
TZa oo. Using Eq. (^Of), we obtain that the free energy of a stripe of length C is 

52 



2an 



(53) 



The term in the brackets characterizes the rigidity of a stripe. As can be seen from Eq. 
this rigidity becomes negative at a critical value of 6 = 6±, where 



5, 



(2(To/t e 



3 xl/2 



(54) 



At 5 > 6± the stripe becomes unstable with respect to wriggling (see Sec. |1VCD 
Taking the limit TZa ^ 00 in Eq. (|5TD , we obtain 



6_ 

en 



(55) 



Note that the stripe solutions of Eq. (|34[) exist only when £s > Ine ^, so to have a solution 
we must have 6 > elne~^ 0, Also, the width of a stripe is limited by Cg ~ e^^, in which 
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case our linearization approximation to Eq. pO| ) breaks down (see also [Q]). According to 
Eq. (|55|), the region of existence of stripes is wider than that of spots. Also, note that for 
5 ^ Sb (when the spot is stable), the width of the stationary stripe Cg ~ g-i/^ ^ This 
deviation from the equilibrium length scale given by Eq. (Ill) is essentially related to the 
one-dimensional nature of the stripe, for which the curvature effects are absent. 

B. Hexagonal and lamellar patterns 

When the spots or stripes are introduced into the system, the basic interaction between 
them is repulsion [see Eq. (pS])]. In an equilibrium configuration, the domains will therefore 
go as far apart from each other as possible. If in the end the distance between them is 
greater than the screening length e~^, essentially they will not interact, so their behavior 
will be that of the solitary patterns discussed in Sec. [IV A| . The situation changes, however, 



when there are so many domains in the system that even in the close-packed arrangement 
the distance between them becomes less than the screening length. This is in fact a generic 
situation that is realized whenever the value of is not close to —1 (or 5 ~ 1). In this case 
the domains will strongly interact with each other, arranging themselves into a multidomain 
pattern, so in order to decrease the energy of the long-range repulsion, the domains not only 
adjust their positions, but also their geometric characteristics. 

Let us consider the simplest of the multidomain patterns in two dimensions, namely, 
the periodic hexagonal and lamellar patterns. The equilibrium characteristics for several 
major types of periodic patterns described by Eq. (^ in the limit e ^ were found by 
Ohta and Kawasaki [|18[. They carried out a rather involved calculation of the free energy 
using the Ewald summation method. Their results can be obtained by the simpler, although 
approximate, Wigner-Seitz method Consider a hexagonal pattern made up of circular 
domains, for example. In such a pattern ip will satisfy Eq. (^) with no flux boundary 
conditions on the boundaries of the hexagonal Wigner-Seitz cell. Instead of solving this 
problem, let us consider a single domain inside a circular cell whose area is equal to the 
area of the Wigner-Seitz cell (a similar approach was used in ||6^). Then Eq. ( p3l) with no 



flux boundary conditions can be easily solved. Furthermore, to the leading order in e the 
screening term e'^n'^ip in Eq. (^) can be neglected, if the period of the pattern Cp ^ e^^. 
This solution can be used to calculate the contribution from the long-range interaction to 
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the free energy by substituting it to Eq. (pTl). Note that the one-dimensional analog of this 
method is exact, so it can also be used to calculate the free energy of the lamellar pattern. 

Let the positive domains in a hexagonal pattern have radius TZg and period Cp < e-\ It 
is convenient to introduce the fraction / of the total area of the system occupied by positive 
domains. The condition of Eq. (|^) implies that this fraction is related to as 

/ ^ (56) 

since inside the domains (away from the interfaces) (p — =tl [see the discussion after Eq. (p5|)]. 
In the hexagonal pattern TZg and Cp are related via 

From this equation one can see that when (p is not close to —1, the values of TZg and Cp 
are comparable. Note that in reality Eq. (0) is approximate, since generally the domains 
forming a hexagonal pattern are not ideally round. However, according to the numerical 
simulations, the deviations from the circular shape are very small when < 0, so one can 
safely assume the domains to be ideally circular all the way up to = 0. 

In Eq. (^) the period of the pattern has not been specified. In fact, an infinite set 
of solutions in the form of hexagonal patterns with different periods exists for — 1 < < 
(asymptotically). All these solutions locally minimize the free energy of the system. 
However, among all hexagonal patterns there is a pattern with a particular period C* for 
which the value of the free energy is the lowest. It is clear that if the asymmetry between 
the positive and negative domains is strong, the domains will tend to form a close-packed 
structure, so for 5 ^ 1 in ci = 2 this pattern is expected to be the global minimizer of the 
free energy. 

Using the Wigner-Seitz method, we find that the period C* of the hexagonal pattern with 
the lowest free energy is (Appendix ^) 

As ought to be expected, T^.^ ~ £* ~ e^^/'^. It is interesting to note that this result agrees 
with the exact calculation of Ohta and Kawasaki in []T8| to within 0.1% for all / < 0.5. 
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Similarly, in the case of the lamellar pattern the period Cp and width Cs of the stripe are 
related as 



2 



(59) 



A calculation analogous to the case of the hexagonal pattern shows that the period of the 
lamellar pattern that has the lowest free energy is given by (Appendix 0) 

1/3 



CI 



-2/3 



60-n 



(60) 



This result agrees with that of |18, 28, 47, |5^, 6^. Let us point out that in one dimension Ren 
and Wei proved that the lamellar patterns of arbitrary periods are the only local minimizers 



of the free energy [47]. 

A comparison of the free energies per unit area of the lowest free energy hexagonal and 



lamellar patterns shows that the lamellar pattern has lower free energy when / > 0.35 [jT8 
For < / < 0.35 the hexagonal pattern has the lowest free energy in two dimensions. 



C. Morphological instabilities of solitary patterns 

An important feature of patterns in systems with competing interactions is the fact 
that under certain conditions they can undergo morphological instabilities which lead to the 
distortions of their shapes and transitions between them 0. In reaction-diffusion systems 



these instabilities have been analyzed in [0, H3, H3, |50|, O, p3|, E7 



Apart from the arguments of Sec. |III B| , the physical reason for the existence of morpho- 



logical instabilities is the fact that the energy of the long-range interaction increases faster 
than the area of the domain as its size gets bigger. Therefore, at some critical size it may 
become energetically favorable for the domain to split into two domains of smaller size or 
significantly change its shape. It is interesting to note that such an instability was first 
analyzed by Lord Rayleigh back in 1882 [|6^] . 

To investigate the morphological instabilities of the domain patterns in systems with long- 
range Coulomb interaction, we start by looking at the simplest possible patterns: spots and 
stripes. This analysis was performed in 0] in the context of reaction-diffusion equations of 
activator-inhibitor type. Here we will re-derive these results using the interfacial approach 
[Eq. 
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1. Spot 



Let us consider a single localized spot first. The fluctuations of the spot's shape are 
the azimuthal distortions of its walls characterized by the azimuthal number m (Fig. 
Because of the radial symmetry, the operator L in Eq. ( ^21) is diagonal in the basis formed 
by the functions e*'"'^, where Lp is the polar angle that represents a point on the interface. 
For a spot of radius TZg ^ e^^ one can neglect screening and use G in the integral for m > 0, 
since these fluctuations do not change the overall area of the spot. Assuming p = e*™"*^ and 
calculating the respective integrals, we obtain 

A = Ao + ^ + ^, (61) 

where Aq is a constant that comes from the curvature and the 2(n- Vip) term. In fact, we do 
not need to calculate this constant from the deflnition. Instead, we can use the translational 
symmetry of the problem and note that A = for m = 1 to flnd that 

Ao = - 2e27^,. (62) 

Note that Aq < and is responsible for the instability of the spot for large enough TZg (see 
also i, H). 

According to Eq. (|61|), a single localized spot becomes unstable (A < 0) with respect to 
the m-th mode when TZg > TZcm, where 

= [-^-^2 -) ■ (63) 

The instability is realized flrst with respect to the fluctuation with m = 2, so the spot is 
always unstable when TZg > 7lc2, where 

7^.2 = (3ao)^/^6-^/l (64) 



The m = case can be treated analogously, this leads once again to Eqs. (^) and (^5|). 
Therefore, comparing Eq. (|6^) with Eq. (^9]), we see that the spot can be stable only when 



6b<6<6c2, (5c2 = 3-^/^KVo^V/3lne-\ (65) 

so the spots are stable only in the limited range of 5 ~ e^^^ <^ 1 (apart from the logarithmic 
terms) . 
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A similar analysis shows that a thin annulus of radius TZa ^ considered in Sec. |IV A 
is always unstable with respect to the m = 2 wriggling mode, so we do not present this 
analysis in detail here. 



2. Stripe 

Let us now turn to the solitary stripe. Let us choose the reference frame in such a way 
that the stripe is oriented along the y-axis in the zy-p\a.ne. Because of the mirror symmetry 
of the stripe in the z-direction, there are two basic types of fluctuations: the symmetric and 
the antisymmetric distortions of the stripe walls, both characterized by the transversal wave 
vector k±_ (Fig. |^). Because of the translational symmetry in the y-direction the operator 
in Eq. ( ^2]) is diagonal in k±. Assuming that p'^ = e*'^-^^ and p~ = ±p+, where are the 
positions of the right and left boundaries of the stripe, respectively, we can calculate the 
integral in Eq. (^2]) at the location of the right wall: 



4e' dS'G,{x-x')p{x') 

^ 26^[l±exp(-/:.v/6V^ + fcl)] ^,,,^ 

where "+" corresponds to the symmetric, and "-" to the antisymmetric fluctuations, re- 
spectively; Cg is the width of the stripe. 

For the stripe, the curvature terms in Eq. (^21) are zero, and the 2{n ■ Vip) term reduces 
to a constant Aq < 0. The case of the symmetric and antisymmetric fluctuations must be 
treated separately. For fc_L ^ e the expression for A = A+ for the symmetric fluctuation is 
(to the leading order in e) 

A, = A„ + .„*i + Hf!l±^5E(zM^. (67) 

K± 

On the other hand, when £s ^ and k±Cs <^ 1, to the leading order in e the expression 
for A = A_ for the antisymmetric fluctuation becomes 



A_ = Ao + (Tokl + 2e^Cs - e^C^^e^^^ + kl. (68) 

Once again, we can use translational symmetry in the z-direction to calculate Aq, since 
A_ = when k± = 0. We get 

Xo = -2e^Cs + e^KCl + 0{e^C^). (69) 
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The analysis of the transcendent Eq. (|67|) with Aq given by the first term in Eq. (|69| 



shows that the instabihty of the stripe with respect to the symmetric distortions of its walls 
(corrugation) with k± = kc occurs at Cg > Cci, where |Q 



k, = 1.13ao-'/'e2/3, = IMa'/'e^'^^ (70) 

According to Eq. (^), the stripe becomes unstable with respect to the antisymmetric 
distortions of its walls (wriggling) at A; — > and Cg > Cc2, where 

C,2 = {2aoKY/h''/'. (71) 

This is also clear from Eq. (|53D. Comparing Eq. (|7T[) and (|70D, one can see that the 
instability with respect to wriggling is realized before the instability with respect to the 
corrugation. In view of Eq. (^5]), the stripe is stable only when 

e\ne-'<5<5^, (72) 

where 5± is defined in Eq. (|5^). Thus, the region of existence of stable stripes is wider than 
that for spots, see Eq. (|65|). 

D. Morphological instabilities of hexagonal and lamellar patterns 



The solution of Eq. (|4^) in the case of an arbitrary multidomain pattern is a formidable 
task. However, a simplification of this problem is possible in the case of periodic patterns. 
Then, by Bloch theorem, the operator L can be partially diagonalized by considering the 
fluctuations modulated by the wave vector k which lies in the first Brillouin zone of the 
underlying domain lattice. The situation here is not unlike the problem of finding the 
band structure of a crystal |^. Below, we consider stability of the hexagonal and lamellar 
patterns in two dimensions. 

1. Hexagonal pattern 

Let us consider a hexagonal pattern of period Cp made of circular domains of radius TZg. 
For each domain centered at R„ let us write the displacement p„ as 

p„(<^) =^a„e*'^-^"-^'"^ (73) 



27 



where the angle represents a point on the interface of each domain. Equation (|^) for the 
fluctuations with a given k in the first Brillouin zone then reduces to 

(~^~ + -^0 - a™ = - ^ Rmm/(k)a„/, (74) 

where m and m' are the azimuthal numbers, Aq is a constant independent of m and k (as- 
suming that with good accuracy n ■ Vip is radially-symmetric in the interface), and Rmm'(k) 
are the k-dependent matrix elements of 4e^Gf^{x — x'). A calculation of Appendix |C 1| shows 
that [H 

167^6^7^, ^ e»(m-mO(^i.+,„+f) 

X J„(|k + k„|7^,) Jw(|k + k„,|7^,), (75) 

where k„ run over the reciprocal lattice, Jm{x) are the Bessel functions of the first kind and 
^k+kn is the angle between the vector k + k„ and the x-axis. 

The value of Aq can be calculated by noting that the translational invariance of the system 
requires that A = for k = and m = 1, so (Appendix |C Ij ) 



Ao = -:S^-2e^7^Jl-^). (76) 

In writing the above equations we assumed that IZs ^ Cp <^ e^^. 

As was shown qualitatively by Kerner and Osipov, for the most dangerous fluctuations 
the wave vector k will lie close to the edge of the Brillouin zone [0]. There are two basic 
types of fluctuations we need to consider: the fluctuations which lead to repumping of the 
order parameter between the neighboring domains (Fig. ^[a]) and the fluctuations which 
lead to the asymmetric distortions of the domains (Fig. ^[b]) ( [0], see also Sec. p.llB| ). 
The analysis of Eq. ( |7^ shows that the most dangerous fluctuations leading to repumping 
have k = — ^2)-, where bi and b2 are the reciprocal lattice vectors which make a 120° 
angle (see Appendix |C 1|) , while the most dangerous fluctuations leading to distortions have 
k = |(bi + b2) [|^. The instability A < occurs with respect to repumping when Cp < Cpo 
or with respect to the asymmetric distortion when Cp > Cp2, where £po,2 depend on e and 
71s/ Cp. The resulting stability diagram obtained by numerical solution of Eq. (|7^ with the 
parameters of the mean-field model from Sec. [II B is presented in Fig. ^. From this figure 
one can see that only the patterns with the period Cpo < Cp < Cp2 are stable. 
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Figure ^ also shows the period of the lowest free energy hexagonal pattern given by Eq. 
(^). One can see that this pattern is stable for all values of IZg/ Cp (except, possibly, for 
IZs/ Cp close to 0.5 where the assumption about the circular shape of the domains ceases to 
be valid). As was noted in Sec. [IV lj| , in two dimensions the lowest free energy hexagonal 
pattern is expected to be the global minimizer of the free energy if / < 0.35, or, equivalently, 
if TZs/ Cp < 0.31, whereas for 0.31 < IZg/ Cp < 0.37 (the second condition means that < 0) 
the global minimizer should be the lamellar pattern. Figure ^, however, does not show the 
transition from the hexagonal to the lamellar pattern, so in fact the lowest energy hexagonal 
pattern is at least metastable for all values of (p at which it exists. 



2. Lamellar pattern 

In the case of the lamellar pattern one should consider the fluctuations that are modulated 
by the wave vector k\\ from the first Brillouin zone in the direction along the normal to the 
stripes and an arbitrary wave vector k± in the transverse direction: 

where are the displacements of the right and left walls of the stripe in the n-th period 
of the lamellar pattern at y = 0. Because of the translational symmetry in the direction 
along the stripes, these fluctuations are the eigenfunction of L in Eq. (^21) . One can then 
reduce the operator L to a 2 x 2 matrix, so after a tedious calculation (Appendix C2, see 
also i H) 

A± = Ao + aokl + R± (fc|| , , (78) 
where R±{k\\,k±) are given by 



R±(fc||, fc±) = — —r-. — X 



/c (l — 2e^'--p COS k\\Cp + e'^^'-p) 
sinh kCp ± [(sinh[A;(£p — Cg)] + cos k\\Cp sinh kCg) 

+ sin^ k\\Cp sinh^ kCg] ^^"^ 



2 



(79) 

where k = a/ e'^K^ + k^. As before, the value of Aq is determined with the aid of the 
translational invariance of the system, which requires that A = for /cy = k± = for the 
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antisymmetric fluctuation. This gives the following value of Aq (to the leading order in e): 

Ao = -26^/:, (l - • (80) 

The analysis of Eq. ([78| ) shows that (in the validity range of Eq. (|79|), that is, when Cp 3> 
In e~^, see also [Q) the repumping instability is not realized for the lamellar pattern. This can 
be explained by a simple argument: the curvature of the stripes is equal to zero, so there is 
no force that would lead to the domain collapse as in the case of the spot. The analysis of Eq. 
([78|) also shows that the most dangerous fluctuation leads to the antisymmetric distortions 



of the stripe and has k\\ = and k± ^ (see also [63, |65[). All other instabilities, such as the 
corrugation instability, occur at higher values of Cp ~ e~^/'^ (compare with [^6|, |6^). Also 
notice that when k± = 0, what corresponds to the one-dimensional situation, the lamellar 
pattern is always stable when Ine^^ <^ Cp <^ e^^ (see Appendix |C 2|) . This is in agreement 



with the result of Ren and Wei that in this situation the lamellar patterns are all local 



minimizers of Eq. (|2^) P? 



Solving Eq. (^) with k\\ = 0, we obtain that the instability is realized when Cp = C*, 
where C* is the period of the lowest free energy lamellar pattern given by Eq. (^) (see 
Appendix |(J 2|) . This result was also obtained by Yeung and Desai in the case / = 0.5 



pq] . Thus, the lowest free energy lamellar pattern is marginally stable with respect to the 
wriggling instability. This fact has a simple geometric interpretation, and should in fact 
be true for any system with long-range competing interactions. Indeed, consider a small 
wriggling modulation of the lamellar pattern (Fig. P). Inside the dashed line rectangle 
the stripes can be considered as straight, but rotated by a small angle i). This pattern can 
be again considered as a lamella, but with a smaller period Cp = Cp cos i). The free energy 
of the system in this case will increase if Cp < C*, since the free energy is a decreasing 
function of Cp for these values of Cp, or decrease if Cp > C*, since there the free energy is 
an increasing function of Cp. The case Cp = C* is marginal. Therefore, the lamellar pattern 
will be unstable with respect to wriggling if Cp > C* or stable otherwise. 

It is interesting to note the analogy between the lamellar patterns and smectic A phases. 
In smectics the long-wave modulations of the layered structure cost free energy AF oc 
(j3k'^^ + Kik^ \uz{k\\,k±)\'^, where Uz{k\\,k±) is the amplitude of the layer displacements 
modulated by wave vecotors k\\ and kj_ along and perpendicular to the layers, respectively 
J^, This is precisely what we get for the pattern with Cp = C* in the limit of small 
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k\\, k± (see Appendix |C 2|) . Furthermore, in view of this analogy the long- wave instability 
of the lamellar patterns with Cp > C* is equivalent Helfrich-Hurault instability of smectics 
under stretching deformations. Also note that under the influence of thermal fluctuations 
the lamellar pattern with Cp = C* is subject to the Landau- Peierls instability We 



would like to point out, however, that all this does not apply to the metastable lamellar 
patterns with Cp < C*, which, according to our calculations, have finite shear modulus. 

V. SCENARIOS OF DOMAIN PATTERN FORMATION 

The analysis in the preceding sections shows that the patterns in systems with long-range 
Coulomb interactions are very sensitive to the parameters and e and can undergo various 
instabilities. As the temperature is lowered, both e and \(f)\ rapidly decrease, see, for example, 
Eqs. (1^) and (^. In this situation a small variation of the temperature may trigger complex 
spatiotemporal behavior in the system. Now we would like to ask the following questions: 
how do the patterns form in the initially homogeneous system and how do the patterns 
already present in the system react to changes in the external parameters? In principle, 
to answer these questions we need to specify the dynamics of patterns. This question is 
quite complicated and significantly depends on a particular system, despite the universality 
possessed by the free energy (for various examples, see [|, |28|, ^). However, if the 
dynamics of the system is dissipative, it will result in the decrease of the free energy of 
the patterns with time. Note that since the free energy functionals in Eqs. (^) or (p9D are 



bounded from below in systems of finite volume, in the absence of the noise the patterns 
must evolve to local minimizers of the free energy. To mimic this behavior, we will use the 
simple gradient descent dynamics defined by 

where F is given by Eq. (|^). This evolution equation is in fact applicable to a number 
of systems with non-conserved order parameter |^, ^, However, our conclusions 



should not qualitatively depend on this particular choice, since the evolution of patterns 
will generally be guided by the free energy landscape and the morphological instabilities of 
the patterns. Note that Eq. ( ^T]) is equivalent to a reaction-diffusion system with the fast 



inhibitor and can be reduced to a free boundary problem in the limit e ^ ||56|, which, in 
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turn, is the gradient descent dynamics for the interfacial free energy 



i5| 



Nucleation 



The first question is how the domain patterns form in the system in the first place. As 
was discussed in Sec. |II at sufficiently high temperatures above Tc the homogeneous phase 



is the only equilibrium state. In the mean-field model of Sec. [II B| the homogeneous phase 
becomes unstable as the temperature of the system is lowered. At = this happens when 
e(r) = ec ~ 1. On the other hand, when the original (unsealed) value of (p is different from 
0, the homogeneous state will remain stable even for lower temperatures. The greater the 
(unsealed) value of |0|, the lower the temperature at which the homogeneous phase will lose 
its stability. This means that for a <^ 1 the homogeneous phase will typically remain stable 
in a range of r for which e <^ 1. 

On the other hand, as was shown in Sec. |1V A| , when the scaled value of |0| is less than 
\(j)b\ = 1 — ^b, in addition to the stable homogeneous phase the system can support stable 
domain patterns (spots). In a narrow range of < 1 01 < \(pb\, where = ^ ~ ^m, the 
spots will be energetically unfavorable. On the other hand, in a wide range of |0| < \(f)m\ the 
domain patterns will have lower free energy than the homogeneous phase. In the mean-field 
model of Sec. [11 B| the homogeneous phase remains stable as long as > |0c| — for 



e ^ 1. Therefore, at |0c| < |0| < |0m| the homogeneous phase is metastable. 

The metastability of the homogeneous phase implies the possibility of nucleation of the 
domain patterns as a result of thermal fluctuations. It is natural to assume that the nucleat- 
ing droplet in two dimensions is a spot (it is localized in space and radially-symmetric). Let 
us consider the nucleation of the positive domains from the negative homogeneous phase. 
In this case the value of is negative, so, as the temperature of the systems decreases, the 
value of (f), as well as 6 [recall Eq. (PD], increases. As was shown in Sec. |1V A| , at 5 > 5^ there 



are two spot solutions. The spot with the radius TZ = TZn is in the unstable equilibrium 
with the homogeneous phase (see Fig. |I|). Therefore, it is this solution that should play the 
role of the nucleation droplet in our system. Note that the radius of the nucleation droplet 
Tin < Tlmin and is bounded for all 6, in contrast to systems with first-order phase transitions. 
This is a distinctive property of the systems with long-range Coulomb interactions. 

According to Eq. (^3]), for 6 close to 6h the free energy cost of the nucleation droplet is 
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(apart from a weak logarithmic dependence) 

AFdrop oc > 1. (82) 

Note that in three dimensions the same arguments give AFdrop oc t'^l"^ ^ 1 also. Since this 
free energy barrier is high, the arguments of the nucleation theory apply here. 

In the narrow range of 5;, < 5 < a nucleation event will result in the formation of a 
stable spot whose free energy is higher than that of the homogeneous phase. Therefore, in 
this situation the spot itself will be metastable and will decay back into the homogeneous 
phase. However, when b is not in the immediate vicinity of 5^, the free energy barrier 
the spot has to overcome to decay will also scale as in Eq. (^), so such spots will be 
long-lived metastable states that can be excited by thermal fluctuations. Therefore, the 
thermodynamic equilibrium state of the system for these parameters is a rarefied gas of 
noninteracting spots. In this situation the spots will play the role of quasipartides. 

When S exceeds 6m, the spot becomes thermodynamically more favorable than the ho- 
mogeneous phase. For 6m < S < 6c2 the spots with the radius TZg are stable, so a nucleation 
event will result in the formation of a single stable spot. This is another distinctive feature 
of nucleation in our system: sufficiently close to 6b a single nucleation event will result in a 
formation of only one spot. However, in order for the system to come to the equilibrium, it 
has to become filled with spots, so the transition from the metastable homogeneous phase 
to the equilibrium multidomain pattern requires many nucleation events. These events will 
occur on an extremely long time scale r^^^ci ~ e'^^-^'^^p 3> Tj-^i, where c > 1 is a constant, Trei 
is a characteristic system's relaxation time. Dynamically, this phenomenon can be identified 
as aging ||70(| . 



According to Eq. p9|), as the temperature decreases, and, therefore, and 6 increase, 
the radius of the nucleation droplet gets smaller (in the scaled units), while the radius of the 
stable spot becomes larger, so at 5 = 6c2 the spot with the radius TZs will become unstable 
with respect to the morphological instability (see Sec. |IV C| ) . In this case the nucleation 
scenario will change. Instead of a single spot, a more complex pattern will form as result of 



a single nucleation event (see Sec. [VB| ). At the same time, the nucleation barrier decreases 
with temperature. At 5 ~ 1 the free energy barrier becomes AFdrop ~ 1 [see Eqs. ( ^31) 
and (^U])]. So, in the renormalized model of Sec. p.1 D| nucleation becomes meaningless for 



these values of 6, and one can talk about the instability of the homogeneous phase. Let us 
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point out that for 5 > 5c2 the annulus should also be considered as a potential candidate for 
the nucleation droplet. The comparison of the free energies of the annulus and the spot of 
radius TZn shows, however, that the spot always has lower free energy. 

Let us consider nucleation in the renormalized model of Sec. |11 D| in more detail. Depend- 
ing on the cooling rate, several situations are possible. When the cooling rate is very small 
(with the characteristic time scale much longer than Tnuci), the system will have enough time 
to equilibrate and will get filled with spots. If the cooling rate is such that its characteristic 
time scale is comparable to Tnuci, the system will enter into the aging regime. Notice that 
these phenomena will occur only in a narrow range of temperatures at which 6 ~ 6b- 

The situation will change qualitatively when 6 > 6c2, when the temperature falls below 
the value at which the spot becomes unstable with respect to a morphological instability. 
In this case a single nucleation event will produce a spot which will further develop into 
a more complex extended pattern, filling up the whole system (see Sec. |VB| ). The time 
scale of this process is Trei and is much shorter compared to Tnuci, so, if the cooling rate is 
sufficiently fast, only a single nucleation event is enough to create a pattern that will fill the 
entire system. This will also be the case in the mean-field model of Sec. |11 B| . Finally, if 
the cooling rate is very fast, the system will not have enough time to nucleate even a single 
domain, so it will enter the region in which the homogeneous state of the system is unstable. 
In that case a pattern consisting of the domains whose characteristic size is comparable 
with the correlation length will form spontaneously and then evolve towards equilibrium via 
coarsening (see Sec. [VQ) . 

If small local inhomogeneities exist in the system, they can work as the nucleation centers. 
One could, for example, have a slightly nonuniform distribution of across the system. If 
their amplitude and size are not very large, the nucleation events will produce stable spots 
that will be pinned to the locations of these inhomogeneities. If, on the other hand, the 
amplitude and size of these inhomogeneities are large enough, the spots that nucleate at 
their locations may be unstable with respect to the morphological instabilities, so they can 
work as the nucleation centers for the spatially extended patterns. 
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B. Growth of complex patterns 



As was discussed in Sec. [V A| , typically a spot that forms as a result of a nucleation event 
will be unstable with respect to the deformations of its shape. After such a spot is formed, 
it will start to grow into a more complex pattern on the time scale of Trei -C Tnuci- Therefore, 
in the process of growth of such a spot thermal fluctuations become unimportant. A typical 
evolution of a pattern in this situation is shown in Fig. |^ which shows the result of the 
numerical solution of Eq. (^Tj) in this parameter regime. At long times, the system evolves 
to a disordered metastable equilibrium pattern. 

The evolution of the unstable localized patterns within the framework of Eq. (^) and the 



corresponding interfacial dynamics problem was studied in detail in |^5|, ^ q3 • The solution 
of the interfacial dynamics equation shows that for sufficiently large 6 the morphological 



instability of a spot will always lead to self-replication of spots [^]. As a result of the 
instability, the spot grows more and more distorted, until at some point the interfaces touch, 
what leads to the pinch-off and splitting (fission) of one domain into two. The daughter 
domains move away from each other, and the process of splitting repeats itself. This self- 
replication process will continue until the whole system gets filled with the multidomain 
pattern O 0]. 

Note that these results hold in the asymptotic limit e — > 0, in which TZg ~ e^"^/^ <^ e^^^ 
so the screening effects may be ignored. On the other hand, for reasonably small but finite 
e the effective interaction may get truncated at distances comparable to the sizes of the 
spots. In particular, as the spots move away from each other after splitting, the interaction 
between the distant portions of their interfaces may get screened, so the spots will remain 
connected by a stripe as they move apart (see Fig. ^ t = 580). For Sc2 < S < 6c3 (the 
latter corresponds to the value of 6 at which the spot becomes unstable with respect to 
the m = 3 mode) the tips on both sides of the stripe will be stable, so as a result of the 
destabilization a spot will transform into a stripe spanning across the system. Note that 



according to Eq. (|53D , this can happen only when 6 > 6±, when the stripe is energetically 
favorable. According to Eqs. ( |54D and (|65D, for not very small values of e it may be possible 
to have 6c2 ~ S±. Note that for these values of 6 the newly grown stripe will destabilize 
with respect to wriggling and fill up the entire space of the system. Furthermore, the stripe 



segments with the highest curvature may become unstable with respect to fingering [57 
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This is what we see in the numerical simulations of Eq. ([81| ) in this parameter range. When 
^ ^ ^c3, the tips of a stripe growing as a result of a splitting event can further destabilize 
with respect to the m = 3 mode, what will result in tip splitting and the formation of a 
labyrinthine pattern. Note that similar results were obtained in the case of reaction-diffusion 



systems with weak activator-inhibitor coupling [HH . 

From the arguments above it is possible to conclude that following a nucleation event at S 
not very far from 6c2 and for not very small values of e the dominating pattern morphology 
is that of the stripe. Nevertheless, for larger values of 6 the spot morphology starts to 
compete with the stripe morphology. This is because for large 6 the interfaces will tend to 



split, so the forming labyrinthine pattern will become disconnected When the value 



of e is decreased, one should find coexistence of both spot and stripe morphologies in the 
patterns formed as a result of the destabilization of a spot (Fig. So, all the processes 
associated with the dynamics of the interfaces: fission, elongation, tip splitting, fingering, 
and wriggling [^, should generally be important for the evolution of a single unstable spot. 

C. Coarsening and disorder 

If the system is quenched deeply into the region in which the homogeneous phase is 
unstable, at first a small-scale multidomain pattern will form. Since the effect of the long- 
range interaction can be seen only on the length scale R > e~^/^, initially the long-range 
interaction will be negligible. Therefore, immediately following the quench the pattern will 
undergo transient coarsening. Note that in systems with the non-conserved order parameter 
this transient coarsening may proceed at arbitrary 0, since the long-range interaction ensures 
the conservation of the total amount of the order parameter. For example, in the case of 
Eq. (|81|) the interfaces of the domains will be driven by curvature subject to global coupling, 
so the characteristic radius of the domains will obey the standard t^^^ law independently 
of the volume fraction |3^, |7T], |7^. As a result, the characteristic size of the domains will 
grow until it becomes comparable with the equilibrium size of Eq. (|T^ . Then the long-range 
interaction will stabilize the pattern, so at some point the coarsening will become arrested 
(see also [^, 73 1). This scenario is observed in the experiments on thin diblock copolymer 



films l?^. Note that this coarsening can be viewed as the consequence of the repumping 



instability discussed in Sees. [II B and IV D 
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When the temperature in the mean-field model from Sec. [II B| is slowly lowered, so that 
the value of |0| becomes lower than |0c|, the homogeneous state becomes unstable with 
respect to the fluctuations with the wave length A ~ e^^^^ (see Sec. [11 B| ). Thus, at the 
threshold of the instability the domain pattern with characteristic size ~ A will start to 
form 0]. These domains will still be smaller than the equilibrium size R ~ e^^^'^, so the 
instability will be followed by coarsening, just like in the case of the renormalized model and 
in the mean-fleld model not close to (pc- 



The results of the simulations of Eqs. (81) displaying transient coarsening are presented 
in Fig. H. In all these simulations the initial conditions were taken as = plus small 
random noise. One can see that the morphology of the pattern is determined by the volume 
fraction of the positive domains. When = the pattern that forms in the end of the 
simulation is a bicontinuous domain patterns similar to patterns forming in the process of 
Ostwald ripening after the critical quench [^. When there is a small asymmetry between 
the positive and negative domains [Fig. |^(b)], the pattern in the end of the simulation looks 
like a collection of disconnected spots and stripes of different shapes and sizes. When the 
asymmetry between the positive and negative domains is strong [Fig. |^(c)], only the spot 
morphology survives, and in the end the pattern is made of a polydisperse mixture of spots. 

Let us emphasize that the patterns that form at the end of the simulations of Fig. |^ 
do not change in time, that is, they are metastable. Each of these patterns is completely 
disordered, and is in no way close to the perfectly ordered patterns that are expected to be 
the global minimizers of the free energy. In the absence of the noise the shape of the pattern 
at long times is determined only by the random initial conditions. The numerical analysis 
of Eqs. (0) shows that by changing the random seed which determines the initial condition 
at the start of the simulation one will get totally different metastable patterns in the end, so 
the system is in fact very sensitive to the initial conditions. This also suggests that besides 
the ordered equilibrium patterns, there exist a huge number of irregular metastable patterns 
which locally minimize the free energy. Thus, a typical pattern that should form as a result 
of the fast quench must be disordered. 

Once the metastable equilibrium is achieved, the patterns will evolve by thermally acti- 
vated processes. Indeed, in order for a pattern with lower free energy to form, some of the 
domains may have to disappear and some may have to be created, since the connectivity of 
different metastable domain pattern is not generally the same. This requires to overcome 
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large free energy barriers AF oc e~^/^. Once again, on the time scale Xnuci ~ e'^'^^ the system 
will enter the aging regime. To mimic this situation, we performed a numerical simulation of 
Eq. (^) with a special initial condition in the form of a metastable hexagonal pattern with 
a single bigger spot in the center (Fig. ||). As time goes on, the pattern tries to adjust to 
accommodate a defect it is presented with. Let us emphasize that, according to Fig. H, the 
defect propagates to distances much larger than the characteristic size of the domains. In 
the end, the pattern becomes completely disordered, with no traces of the original hexagonal 
ordering. 

Let us now ask a different question. Suppose that we already have an equilibrium domain 
pattern in the system. What happens if at some moment the temperature of the system is 
raised or lowered? This question is related to what happens if the system is gradually cooled 
below the transition temperature. Suppose the system is initially occupied by the lowest 
energy hexagonal pattern (at a given temperature). When the temperature of the system 
is lowered, the values of e and \(f)\ will decrease, so this pattern will no longer correspond to 
the equilibrium pattern. To see this, let us write down the length scales in Eq. (|l^) in the 



original (unsealed) variables. Using the mean- field scaling |42| and the definition of e from 
Eq. (H), for example, we obtain 

L~a-V2|r|i/2. (83) 

One can see from here that if the temperature is abruptly lowered, the equilibrium size R 
will increase. At the same time, the physical size of the pattern will remain the same, so the 
relative size of the pattern will decrease with respect to the new value of R. 

If the temperature drop is sufficiently small, the pattern will remain metastable (see 
Sec. [IV D| ). However, when the temperature falls below a certain critical temperature, the 
pattern will become unstable with respect to repumping (Sec. [IV D|) . The repumping will 



lead to the collapse of a fraction of the domains and growth of the rest, so effectively, this 
will be equivalent to the increase of the characteristic interdomain distance. The resulting 
pattern will again be metastable. Note however, that it will necessarily become disordered, 
since randomness is involved in the destabilization of the hexagonal pattern. When the 
temperature gets lower, the metastable pattern will again destabilize, and produce a new 
metastable pattern with a greater characteristic domain size. This process will go on. Thus, 
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we will have a stepwise relaxation process creating disordered patterns. Similar effect will 
be realized if one takes a hexagonal pattern as an initial condition and gradually raises 
the temperature. At some moment the pattern will become unstable with respect to the 
asymmetric deformations, so domains of complex shapes will start to form, thus effectively 
making the distance between the patterns smaller. Such a metastable domain pattern will 
further destabilize at higher temperature. These conclusions are confirmed by the numerical 



simulations of Eqs. (|8l|) [Q. Note that these arguments imply that disordered patterns 
will form even as a result of a slow (but fast compared to t^uci) quench below the transition 
temperature. All this indicates that disorder is an intrinsic state of the systems with long- 
range competing interactions. This is also seen in the experiments 75, 76 . 

Similarly, when one starts with a lamellar pattern and raises the temperature, the pattern 



will become unstable with respect to wriggling (see Sec. IV D ). If the temperature is further 
increased, the corrugation instability and fingering will follow. Notice that in contrast 
to the hexagonal patterns, the lamellar pattern will always remain metastable when the 
temperature is lowered, since there is no repumping instability in this case. This means that 
a metastable lamellar pattern is more likely to survive after a slow (but fast compared with 
Tnuci) critical quench. 



VI. CONCLUSIONS 



In this work we have presented an energetic approach to the study of inhomogeneous 
states (patterns) in systems with competing short-range attractive and long-range repulsive 
Coulomb interactions. Our approach becomes universal for systems with weak Coulomb 
interaction in the vicinity of the microphase separation transition, thus allowing to treat a 
variety of physical situations which involve competing Coulomb interactions. 

By the very definition of the domain patterns, the width of the domain wall should be 
much smaller than the characteristic size of an individual domain. This requires that the 
Coulomb interaction be sufficiently weak in order for these patterns to be feasible. On the 
other hand, one can take advantage of this and study these patterns in the asymptotic 
limit of infinitely weak Coulomb interaction. This poses a challenge, however, since this 
interaction is a singular perturbation to the short-range interaction. 

We have performed an asymptotic analysis of the free energy in the limit of vanishingly 
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small strength of Coulomb interaction (e — > 0). Our main finding is that in this limit the 
energetics of the patterns are described by the locations of the domain interfaces. In fact, 
an important hierarchy of length scales appears in the system [Eq. (|16])]. Our second major 
observation is that the characteristic size of the domains in a stable domain pattern has to 
scale as e~^/'^. This is different from the similar estimates based on the properties of global 
minimizers of the free energy []TB|, BR BR What we showed in Sec. [llllj| in general and 



Sees. [IV Q and |1VD| for particular patterns is that unless this scaling is obeyed, the pattern 
cannot be a local minimizer and thus thermodynamically stable. 

In our analysis, the starting point was the mean-field free energy functional from Eq. (P. 
We chose to perform our calculations using Eq. (j^) for two reasons. First, this is a univer- 
sal functional that is obtained in the vicinity of the microphase separation transition and 
therefore may be applied to a variety of systems. Second, using this functional we could 
obtain very explicit results, making our presentation more tractable. It is not difficult to 
see that all our calculations can be extended to the more general functional from Eq. (p. 
The only difference is that in the case of Eq. (||) the "positive" and "negative" domains may 
have asymmetric linear response coefficients k± instead of a single k in the case of Eq. (^. 
Nevertheless, in the case of Eq. ([1|) we can choose k = n^f + — /), where / is the 
volume fraction of the "positive" domains (similar ideas were used in [^). Indeed, since 
K, is responsible for screening, we can average the response of the order parameter on the 
scale of the domains, which is much smaller than the screening length (Sec. |11 lj|) . The new 
definition of k, also takes into account that to the leading order the (locally) averaged value 
of the volume fraction / is independent of space. The latter can be easily seen from the 
analog of Eq. (|36D obtained from Eq. (|l]), if one integrates this equation over a closed volume 
of size e~^/^ ^ ^ ^ e~^, uses the Gauss theorem, and takes into account that (p is nearly 
constant in the domains and iVipl ^ ip/R ^ e^/^ (see Sec. [II C|) . 

Let us point out that the interfacial representation of the free energy given by Eq. (p9[ ) 
which we obtained from the free energy functional in Eq. (^ in the asymptotic limit e — > 
may in fact itself form a basis for studying the domain patterns in systems with long-range 
Coulomb interactions, provided that the driving force for the formation of these patterns 
is the competition of the Coulomb energy with the surface energy (see, for example, 
r3[ , [55[]). In this formulation our results can be applied to an even wider range of systems. 



which may not generally possess a free energy functional, like the one in Eq. (^. For 
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example, our asymptotic results should apply to the ferromagnetic nearest-neighbor Ising 
models frustrated by Coulomb interactions ^ ^ ^ We argued that in these 
systems thermal fluctuations should only renormalize the effective coupling constant of the 
Coulomb interactions without qualitatively affecting the overall picture (Sec. [11 L)| ). These 
predictions are difficult to compare with the recent Monte-Carlo simulations |S^, 0] because 
of the limitation of the latter on the system's size. Nevertheless, the result of P3] about the 
location of the microphase separation transition, which gives Tc ~ a^, with x ~ 0.25 — 0.35, 
is not far from our prediction from Sec. II D of ~ a^-^^ for the three-dimensional Ising 
model. Note that we do not expect to find an avoided critical behavior discussed in the 



context of the mean-spherical models [77 



An interesting question arising in systems with long-range competing interactions is the 
nature of the thermodynamic ground state below the microphase separation transition tem- 
perature. We emphasize that our stability analysis of stationary patterns only addresses 
small-scale thermal fluctuations, so we are really talking about metastahility of these pat- 
terns. At the same time, rare large-scale thermal fluctuations may lead to nucleation or 
transitions between different metastable patterns (see Sees. V A and |V C| ). In this sense, 
if there are enough metastable patterns, the global minimizer of the interfacial free energy, 
which is presumably a highly symmetric periodic pattern |TB[, has little to do with the 
thermodynamic ground state of the system. 

In fact, we see that the stationary metastable patterns that form in one way or another 
are typically highly disordered (Sec. [V Q) . Although the basic interaction between different 
domains is repulsion, the domains rarely arrange themselves in a close-packed fashion. The 
reason for that is that even though the interaction between the domains is repulsive, the 
range of this interaction, which is determined by the screening length L ~ is much 
greater than the characteristic interdomain distance R ~ e~^/^. So, a single domain interacts 
simultaneously with many other domains and not only with its nearest neighbors. Therefore, 
the optimization of the free energy becomes a collective problem, and a huge number of 
disordered metastable states appears. Then, the configurational entropy of these metastable 
disordered states may overwhelm their energy disadvantage [[70| , [7^ . 

Furthermore, a large-scale fluctuation whose size is comparable with the domain size 
~ R may propagate its action to the much larger distance ~ L (see Fig. It would 
seem natural to expect that even if the system is in the lowest energy state, a sufficiently 
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strong fluctuation will frustrate a region much larger than the size of such a fluctuation, 
which may lead to increase of the degree of disorder with time. In this sense systems with 
long-range competing Coulomb interactions can be considered as structural glasses pO, B? 



We emphasize that in these systems the disorder is self-induced. As we showed in Sec. |V (J 
these systems can age on very long time scales and exhibit complex relaxation phenomena 
even in the case when the equations of motion for the patterns are very simple. Note that 
in a recent paper, Schmalian and Wolynes came to similar conclusions on the basis of their 



replica analysis of Eq. treated as an effective hamiltonian [^. Their calculations suggest 
that the number of metastable states grows exponentially with the system's volume, leading 
to an ideal glass transition below the microphase separation transition temperature. Also 
note that spin systems frustrated by Coulomb interactions were proposed for studying glassy 



behavior in the supercooled liquids ||79| 
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APPENDIX A: FREE ENERGY 

Here we present the details of our manipulations of the interfacial free energy from 
Eq. (^) and the effective field ip from Eq. (|2^) . Let us first show the derivation of Eq. (^Sf) 
from Eq. (^). In view of Eq. (p^ , we have 

1 r,. 



-^long-range = 2 J ^ ^^^""^ " '^^'^ ' ^^^^ 

According to Eq. with 0sh = il? "we have 

+2t^ [ d'^x'G.ix ~ x'), (A2) 



n4 



where the first integral is over the whole space. This integral is equal to l/(eK)^, according 
to Eq. ( p5D . Substituting this back into Eq. ( |A1|) , after simple algebra we arrive at Eq. ( p8D . 
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Let us now derive Eq. (pQ]) from Eq. (pSD- Using Eq. p5| ) and applying Gauss's theorem, 
we calculate the long-range term: 

[ f d'^xd'^x'G,{x-x') = 

2 



[ [ d'^xd'^x' {5^'^\x-x') + V^G,ix-x')) 

4 / d^x + ^ f rfV I dS{h ■ VG,{x - x')} 
Jn+ 1^ Jq+ J 



2 

K^d 



dS{fi ■ x) 
(fdS [ d'^x'V ■ {nG,{x - x')} = 

— r— (t dSifi ■ x) (D dS (D dS'ih ■ n)GAx — x'). 

K^d J J J 



(A3) 



Let us now derive Eq. ( p9D . Using Eq. ( p5D , Eq. ( [A2| ) , and Eq. to express the (5-function 
in terms of G, we get 

1 + 



+4/ d'^x'V^{GAx-x')-G{x-x')} = 

+ -2 idS'{h' ■ V'(a - G)}, (A4) 
n 1^ J 

where we applied the Gauss's theorem. 

Now let us calculate the first and second variations of the interfacial free energy. Let 
p{x) be a normal displacement of the interface at point x on the interface, which is positive 
if the displacement is into the positive domain and vice versa. Note that according to our 
definition, p > corresponds to shrinking of 

Up to second order in p, the change of the surface free energy AFgurf is given by a 
well-known formula (see, for example. 



AFsurf = -2(To i dSHp + 



2 



^dS{\\/±p\^ + 2Kp^), (A5) 



where H and K are mean and Gaussian curvatures at a given point of the interface, re- 
spectively, and denotes the gradient restricted to the interface. The mean curvature 
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is positive if the positive domain is convex. The change of the long-range contributions to 
the free energy is given by an integral over a thin layer of thickness p over the interface. 
According to Eq. (pH]), we have 



2(1 + 



^-^long— range 

dS / dz{l-2Hz + Kz^) 
Jo 

p{x) 



dS J dz{l - 2Hz + Kz^) x 
d'^x' GJx - hz - x' 



+ 



+2e^^dsj^^ dz{l-2H{x)z + K{x)z^)x 
j dS' j^^^ ^ dz'{l - 2H{x')z' + K{x')z''^) X 



Ge{x — hz — x' + n z')^ (A6) 

where we used the fact that with our definition of the sign of principal curvatures d'^x = 
{1 — kiz){l ~ k2z)dzdS = {l — 2Hz + Kz^)dzdS at a point distance z away from the interface. 
Retaining only the terms up to second order in p and using Eq. (|A2|) , we obtain 

-2 i i 



AFiong-range = -2 j> dS Ipp 

+2 j dS Hi)p^ + jdS{h- VV^)p2 



+2t^ jdS j dS'G, {x - x')p{x)p{x'), ( A7) 

where we expanded ip in the Taylor series in z. From this, and Eq. (|A5|), we get 

dF = -2^ dS{aoH + ^)p, (A8) 
so the critical points must satisfy Eq. (^8|). Similarly, using Eq. (|38| ) in Eq. (|A7|) , we obtain 

S^F = (To jdS (I V±p|^ + 2Kp^) 
+ j dS{2{fi ■ W) - 4(To^^}p^ 

which in view of Eq. (|39|) coincides with Eq. (|41| ) . 
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APPENDIX B: OPTIMAL PERIOD OF HEXAGONAL AND LAMELLAR PAT- 
TERNS 



Here, we give the Wigner-Seitz calculation of the period of the hexagonal and lamellar 
patterns. 



1. Hexagonal pattern 



We start with a hexagonal pattern. Consider Eq. (|23D on a disk of radius TZ = 
3^/^£p/-\/27r, with no flux boundary conditions. Neglecting the term e^K^ip and usin^ 
Eq. (^61), we write 



+ 



+ e'{e{n, - r) - e{r - 7^,) - 2/ + 1} = 0, 



(Bl) 



dr'^ r dr 

where r is the radial coordinate and 9{x) is the Heaviside step. The solution of this equation 
that satisfies Eq. (|38|) is given by 



'e^if - ly + - /)} - ko7^-7-^/^ 



< r < 7^„ 



(B2) 



lg2|j^2 _ p^2 ^ fn^ii^f-Ti^ - 21nr)} - ^aoU'^f-^/^ ns<r<n. 



where we took into account that / = IZl/lZ'^. According to Eq. (P?]) with 0sh 
Eq. (^), the long-range contribution to the free energy can be computed as 



±1 and 



-'^ long— range 

(1 - /) / ijd'x -f [ ^d'x = 
Jn+ Jfi- 

27r(l — f) ril){r)dr — 271 f / rip(r)dr. 

Jo JTls 

Combining this with the surface energy Fgurf = ^TxaolZg and using Eq. 
free energy per unit area is 

F 2ao/i/2 e^n^p 



+ 



■(/-I -In/), 



(B3) 

, we get that the 
(B4) 



7r7^2 7^ ■ 2 

so, minimizing this expression with respect to TZ with fixed /, we obtain that the minimum 
is attained aXlZ = 1Z* , where 



2/3 J-1/2 



2(To 



/-ln/-l 



1/3 



(B5) 
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Using the definition of TZ in terms of Cp, this equation is rewritten as Eq. 
2. Lamellar pattern 

Similarly, for the lamellar pattern centered at a; = we get 



e\f - l)x' + \e'f{l - f)Cl 0<x< 



il) = < 



(B6) 



e'fx' - e'fCpX + \e^p{2 - f)Cl CJ2 <x< Cp/2, 
where we used the fact that / = Lsj Hp. Calculating the free energy per unit length, we get 

£ = ^ + ^^^^Jf(l-/)^. (B7) 
Minimizing this expression with respect to £p, we obtain Eq. 



APPENDIX C: STABILITY OF HEXAGONAL AND LAMELLAR PATTERNS 

Here we present the details of our calculations of Eqs. (0) and (|78D. 

1. Hexagonal pattern 

We begin with the hexagonal pattern. We define 

Rmm' (k) = / I d^' 

„ Jo Jo 

^g*mc-imV'+ik.R„^^(j^^ + r(<^') - r(v9)), (CI) 

where = [TZg cos ip, TZg sin ip) and the summation is over the lattice: = niai + n2a2, 
where ai = |£p(\/3, 1) and a2 = |£p(v^, —1). Using the Fourier representation of G^{x — 
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x'), we obtain 

"277 r2TT 



Rmm' (k) = > / "V^ / 

71" „ Jo Jo 



dk.' e«k-R„-ik'-(R„+r(</3')-i"(<^')) 

(27r)2 |k'|2 + e2fi:2 



(C2) 



n 



.„ JO JO 

(C3) 



gi(k+k„).(r{v.)-r(v.')) 



where v = Cp\^/2 is the area of the Wigner-Seitz cell, the sum is now over the reciprocal 
lattice, and we took into account that J2n e*^''"''') -'^" = {Att'^/v) J2n ^i^' — k — k„). 
We proceed: 



|k + k„|2 + e2«;2 



r\mm'(J^) — / , ■- - „ „ X 

TTV 

n 

^^gim^+i(k+k„).r(v) / ^^/g 
JO 



E 



nv Ik + k„ |2 + e2K2 

27r 

^^^irnip+i\l<.+lc„\TZs cos(<^-'i?k+k„ ) 

2tt 







^^/g-jmV-i|k+k„|-R,s cos(<^'-i9k+k„) 
V ^ |k + k„|2 + £2^2 

J^(|k + k„|7^,)^-'"' J™.(|k + kJ7^,), 



(C4) 



where we introduced the angle "i^k+k,, between the vector k + k„ and the x-axis and used the 
integral representation of the Bessel function. After a few algebraic manipulations, this equa- 
tion can be converted to Eq. ([f5|). Using the reciprocal lattice vectors bi = 27r£~^(l/-\/3, 1) 
and b2 = 27r£p ^(l/y^, —1), so k„ = nibi + n2b2, this sums can be evaluated numerically 
by truncating the summation at sufficiently large |ni| and \n2\- 

An alternative representation for Rmm'(k), which allows to explicitly calculate its diagonal 
elements, can be obtained by performing the summation in real space, rather than over the 
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reciprocal lattice. We rewrite Eq. ( P2|) as 

qdq 



Jo 



X 



^f;^g-ig|R„|cos(,9-i?„) ^ 

27r 

27r T 

o A 

^V^' ^-im'if'-iqlZs cos(i/3'-i9) 











/■oo 
■''0 



qdq 

X 







g2 + e^K^ 

/t 

27r 7 q 

"''^^-j<j|R„|cos(i?-i?„)+j(m~m')i? ^ 

^"V„(g7^,)J„,(g7^,), (C5) 



where is the angle between R„ and the x-axis. Calculating the integral over {}, we obtain 

n 

f°° qdq 

/ ~, —Jm-m'{q\'R'n\)Jm{q'R's)Jm'{q'R-s)- (C6) 

Jo q + e K 

In calculating Rmm'(k), to the leading order in e one can neglect the term e'^n^ in the 



denominator of Eq. (^) or Eq. (|C6|). Setting e to zero, we can calculate the diagonal 
elements Rmm(k) for m > 2. After some algebra 

RmM = (C7) 

m 



where we took into account that the integrals in Eq. ( |(J6| ) all vanish for R,„ 7^ 0. Thus, the 
diagonal elements Rmm(k) are independent of k and coincide with those of a single spot. 

Caution, however, is necessary when |k| < e. In this case the k„ = contribution to the 
sum in Eq. (|75D will be singular for m, m' = 0, ±1. Taking only the contribution of k„ = 0, 
for |k| <^ 1 we obtain 



£2^3 \\k\^ + e^^\ 

where we expanded the Bessel functions in the Taylor series and retained only the leading 



term. Now, to calculate Ri,i(k) for k = 0, note that if one formally sets e = in Eq. (75) 



with m = m' = 1, one should get the result of Eq. (P7|). On the other hand, for k = 
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the term with k„ = does not contribute, while for other k„ the term e^K^ in Eq. (^) is a 
regular perturbation and can be neglected. So to calculate Ri_i(k) in the limit |k| 0, one 
has to subtract the e ^ limit of the k„ = term in Eq. ([75| ) from Eq. ( |C7| ). As a result, 
we obtain 

Rii(0) = 2e'7e,-^^-^, (C9) 

where we expanded the Bessel functions in the Taylor series and retained only the leading 
term. 

Since Ro,o(k) 2> |Rmm'(k)| for m, m' 7^ and small |k|, the m = mode is the eigen- 
function of the operator L in Eq. (^2]) for vanishing |k|. The analysis of Eq. (|74D with 
m = m' = and Ro,o(k) from Eq. ( |C8| ) shows that the hexagonal pattern is stable with 
respect to the long-wave modulation of the spots' radii, as long as Cp is large enough. 



2. Lamellar pattern 

Let us now turn to the lamellar pattern. Calculate the matrix elements of 4e^Ge(x — x') 
between the right ( "+" ) and left ( "-" ) walls of the stripe in the zeroth period for a given 
modulation: 



2^2 1^- 

T 

n=— 00 



and 



{+\Ae^G,\+) = {-\Ae^G,\-) 
+00 

^ ^ g— fc|£pn|+ifc||7i£p 

a=— 00 
2g2gifc||£p (^^2kCp _ 



/+|4e2G,|-) = (-|4e2G,|+)* 

2 2 +00 

_ \ ^ —k\Cpn—Cs\+ik»nCp 

~ k ^ 

n=—oo 

+ 



gA:£p _ gifcy^p ^{k+ik^\)Cp — \ I ' 



(CIO) 



(Cll) 



where we introduced k = -^Z e^K^ + /c^, used the fact that the Fourier transform of the 
Green's function in the transverse direction is given by exp[—k\z — z'\)/2k, and summed 
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geometric series. After some algebra, the 2x2 matrix R{k\\,k±) formed by these matrix 
elements can be transformed into the following form 



"' '^^ A; (l — 2e'^'^p COS A;||£p + e^^^p) 

sinh k£p sinh[fc(£p — Cg)] + e*''"^" sinh kCg 

sinh[A;(£p - + 6''''^^^ sinh kC^ sinh kCp 

This matrix can be easily diagonalized, after a few manipulations we arrive at Eq. ([79|). 
Then, Eq. (|80|) is obtained by setting Aq = — R_(0,0) and taking only the leading order 
terms. Note that for k\\ = or k\\ = tt/ Cp the fluctuations corresponding to X± are the 
symmetric and antisymmetric deformations of stripes. 

To obtain the energy of the long-wave distortions of the lamellar pattern, we expand 
Eq. (ffSD with R_ into a series in fcy and kj_ and retain the terms up to quadratic in k\\ and 
forth order in kj_. Then, to the leading order in e, we get 

X-^\f\i-f)Wpkl 

+ (^0 - If {I - ffe'ci^ kl 
+ 3^ [/'(I - + 2/ - 2f )6^/:a fcl, (C13) 

where we used Cg/ Cp = f. One can see that at Cp = C* given by Eq. (pO| ) the coefficient of 
k'j_ changes sign from positive at Cp < C* to negative at Cp > C*, signifying an instability. 
At the same time, the coefficient of k\ is positive for all < / < 1. 

Let us now discuss the stability of the lamellar patterns in one dimension, which can be 
studied by looking at Eq. (|7^) with k± = 0. Setting k± = and expanding in e, after some 
algebra we obtain that to the leading order 

+ f + /(!-/) cos fc|,/:p 
-y/f- 2/(1 COS . 

(C14) 

It is not difficult to verify that according to this equation A_ > for all values of k\\, so the 
lamellar pattern is always stable regardless of the modulation vector /cy in one dimension. 
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A_ 



2e'C^ 



cos knCp 



This conclusion is applicable when lne~^ <^ Cp <^ e~^, when the assumptions of the above 
equations are valid. Note that for Cp outside this range the one dimensional lamellar patterns 
(strata) may undergo a number of instabilities 0. 



[1] M. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 851 (1993). 

[2] B. S. Kerner and V. V. Osipov, Autosolitons: a New Approach to Problems of Self- 
Organization and Turbulence (Kluwer, Dordrecht, 1994). 

[3] A. S. Mikhailov, Foundations of Synergetics (Springer- Verlag, Berlin, 1990). 

[4] R. Kapral and K. Showalter, eds.. Chemical waves and patterns (Kluwer, Dordrecht, 1995). 

[5] F. J. Niedernostheide, ed.. Nonlinear dynamics and pattern formation in semiconductors and 
devices (Springer, Berlin, 1994). 

[6] J. D. Murray, Mathematical Biology (Springer- Verlag, Berlin, 1989). 

[7] M. Seul and D. Andelman, Science 267, 476 (1995). 

[8] T. H. O'Dell, Ferromagnetodynamics: the dynamics of magnetic bubbles, domains, and domain 

walls (Wiley New York, 1981). 
[9] N. Setter and E. L. Colla, eds.. Ferroelectric ceramics: tutorial reviews, theory, processing and 
applications (Birkhauser, Basel, 1993). 
[10] F. S. Bates and G. H. Fredrickson, Annu. Rev. Phys. Chem. 41, 525 (1990). 
[11] A. Y. Grosberg and A. R. Khokhlov, Statistical Physics of Macromolecules (AIP Press, New 
York, 1994). 

[12] F. S. Bates and G. H. Fredrickson, Physics Today 52, 32 (1999). 
[13] C. M. Care and N. H. March, Adv. Phys. 24, 101 (1975). 
[14] V. D. Kovalev and E. L. Nagaev, Sov. Phys. Uspekhi 29, 297 (1986). 
[15] E. L. Nagaev, Phys. Uspekhi 38, 497 (1995). 

[16] R. P. Huebener, Magnetic flux structures in superconductors (Springer- Verlag, Berlin, 1979). 

[17] L. Leibler, Macromolecules 13, 1602 (1980). 

[18] T. Ohta and K. Kawasaki, Macromolecules 19, 2621 (1986). 

[19] F. Liu and N. Goldenfeld, Phys. Rev. A 39, 4805 (1989). 

[20] I. Y. Erukhimovich and A. R. Khokhlov, Polymer Sci. 35, 1522 (1993). 

[21] V. Y. Borue and I. Y. Erukhimovich, Macromolecules 21, 3240 (1988). 

51 



[22] I. A. Nyrkova, A. R. Khokhlov, and M. Doi, Macromolecules 27, 4220 (1994). 
[23] P. G. de Gennes, J. de Physique - Lett. 40, 69 (1979). 
[24] F. H. Stillinger, J. Chem. Phys. 78, 4654 (1983). 

[25] L. Q. Chen and A. G. Khachaturyan, Phys. Rev. Lett. 70, 1477 (1993). 

[26] S. Glotzer, E. A. Di Marzio, and M. Muthukumar, Phys. Rev. Lett. 74, 2034 (1995). 

[27] R. F. Mamin, JETP Lett. 60, 52 (1994). 

[28] C. Yeung and R. C. Desai, Phys. Rev. E 49, 2096 (1994). 

[29] M. Hildebrand, M. Kuperman, H. Wio, A. S. Mikhailov, and G. Ertl, Phys. Rev. Lett. 83, 

1475 (1999). 

[30] V. J. Emery and S. A. Kivelson, Physica C 209, 597 (1993). 

[31] S. Lundqvist and N. H. March, eds.. Theory of inhomogeneous electron gas (Plenum Press, 
New York, 1983). 

[32] U. Low, V. J. Emery, K. Fabricius, and S. A. Kivelson, Phys. Rev. Lett. 72, 1918 (1994). 

[33] P. Viot and G. Tarjus, Europhys. Lett. 44, 423 (1998). 

[34] M. Grousson, G. Tarjus, and P. Viot, Phys. Rev. E 64, 036109 (2001). 

[35] J. Lorenzana, C. Castellani, and C. Di Castro, Phys. Rev. B 64, 235127 (2001). 

[36] J. SchmaUan and P. G. Wolynes, Phys. Rev. Lett. 85, 836 (2000). 

[37] C. B. Muratov, Theory of domain patterns in systems with long-range interactions of Coulom- 

bic type (Ph. D. Thesis, Boston University, 1998). 
[38] W. L. McMillian, Phys. Rev. B 12, 1187 (1975). 
[39] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 
[40] J. R. Smith, Phys. Rev. 181, 522 (1969). 
[41] A. J. Bray, Adv. Phys. 43, 357 (1994). 

[42] L. D. Landau and E. M. Lifshits, Course of Theoretical Physics, vol. 5 (Pergamon Press, 

London, 1980). 
[43] R. Choksi, J. Nonlinear Sci. 11, 223 (2001). 

[44] D. M. Petrich and R. E. Goldstein, Phys. Rev. Lett. 72, 1120 (1994). 
[45] R. E. Goldstein, D. J. Muraki, and D. M. Petrich, Phys. Rev. E 53, 3933 (1996). 
[46] L. D. Landau and E. M. Lifshits, Course of Theoretical Physics, vol. 8 (Pergamon Press, 
London, 1984). 

[47] X. F. Ren and J. C. Wei, SIAM J. Math. Anal. 31, 909 (2000). 

52 



[48] C. B. Muratov, Phys. Rev. Lett. 78, 3149 (1997). 
[49] T. Ohta, A. Ito, and A. Tetsuka, Phys. Rev. A 42, 3225 (1990). 
[50] C. B. Muratov and V. V. Osipov, Phys. Rev. E 53, 3101 (1996). 
[51] T. Ohta, M. Mimura, and R. Kobayashi, Physica D 34, 115 (1989). 
[52] Y. Nishiura and H. Suzuki, SIAM J. Appl. Math. 29, 1087 (1998). 
[53] A. Hagberg and E. Meron, Phys. Rev. Lett. 72, 2494 (1994). 

[54] A. B. Vasileva, V. F. Butuzov, and L. V. Kalachev, The Boundary Function Method for Sin- 
gular Perturbation Problems (Society for Industrial and Applied Mathematics, Philadelphia, 
1995). 

[55] P. C. Fife, Dynamics of Internal Layers and Diffusive Interfaces (Society for Industrial and 

Applied Mathematics, Philadelphia, 1988). 
[56] C. B. Muratov, Phys. Rev. E 54, 3369 (1996). 
[57] C. B. Muratov and V. V. Osipov, Phys. Rev. E 54, 4860 (1996). 
[58] I. Onishi, Y. Nishiura, M. Imai, and Y. Matsushita, Chaos 9, 329 (1999). 
[59] Y. Nishiura and I. Ohnishi, Physica D 84, 31 (1995). 
[60] C. B. Muratov, Phys. Rev. E 55, 1463 (1997). 

[61] J. M. Ziman, Principles of the Theory of Solids (Cambridge University Press, Cambridge, 
1972). 

[62] S. Miiller, Calc. Var. Part. Dif. 1, 169 (1993). 
[63] L. M. Pismen, J. Chem. Phys. 101, 3135 (1994). 
[64] Lord Rayleigh, Phil. Mag. 14, 184 (1882). 
[65] X. F. Ren and J. C. Wei (unpublished). 

[66] S. Qi and Z.-G. Wang, Macromoleculcs 30, 4491 (1997). 

[67] C. Yeung, A.-C. Shi, J. Noolandi, and R. C. Desai, Macromol. Theory Simul. 5, 291 (1996). 
[68] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Oxford University Press, 
Oxford, 1993). 

[69] Y. Oono and M. Bahiana, Phys. Rev. Lett. 61, 1109 (1988). 

[70] J. -P. Bouchaud, L. F. Cugliandolo, J. Kurchan, and M. Mezard, in Spin glasses and random 
fields, edited by A. P. Young (World Scientific, Singapore, 1998), available at LANL archive: 
cond-mat/9702070. 

[71] B. Meerson and P. V. Sasorov, Phys. Rev. E 53, 3491 (1996). 

53 



[72] J. Rubinstein and P. Sternberg, IMA J. Appl. Math. 48, 249 (1992). 
[73] M. Bahiana and Y. Oono, Phys. Rev. A 41, 6763 (1990). 
[74] G. Coulon, B. Collin, and D. Chatenay, J. de Physique II 3, 697 (1993). 
[75] J. A. Cape and G. W. Lehman, J. Appl. Phys. 42, 5732 (1971). 
[76] W. M. Heckl and H. Mohwald, Ber. Bunsenges. Phys. Chem. 90, 1159 (1986). 
[77] L. Chayes, V. J. Emery, S. A. Kivelson, Z. Nussinov, and G. Tarjus, Physica A 225, 129 
(1996). 

[78] D. Chowdhury, Spin glasses and other frustrated systems (Princeton University Press, Prince- 
ton, 1986). 

[79] D. Kivelson, S. A. Kivelson, X. Zhao, Z. Nussinov, and G. Tarjus, Physica A 219, 27 (1995). 
[80] M. E. Taylor, Partial Differential Equations III: Nonlinear Equations (Springer- Verlag, Berlin, 
1996). 



FIGURES 




FIG. 1: The free energy of a spot for different values of (j). The plot of F{JZ) from Eq. ([43|) 
with e = 0.001, ctq = 2\/2/3, and k = l/\/2- For these values of the parameters 5i, = 0.0088 and 
5m = 0.0186. 
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FIG. 2: The free energy of an annulus as a function of TZ for C = Ca given by Eq. (plD, obtained 
from Eq. (|o|) with e = 0.001, uq = 2^2/3, and k = 1/V2. 



m = 2 



corrugation 



wriggling 




m = 3 




FIG. 3: Morphological instabilities of spots and stripes, 
a) 
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asymmetric distortion 







FIG. 4: Two major types of instabilities of the hexagonal pattern. 
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FIG. 5: The stability diagram for the hexagonal pattern with ctq = 2-v/2/3. The dashed line 
corresponds to the hexagonal pattern with the lowest free energy given by Eq. (|58|). 



FIG. 6: Wriggled lamellar pattern. 
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FIG. 7: Formation of a complex pattern as a result of an instability of a spot. Results of the 
numerical solution of Eqs. (|8l|) with e = 0.025 and = —0.6, with no-flux boundary conditions. 
The system is 400 x 460. 
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FIG. 8: Coarsening of the domain patterns at different values of 0: = (a), (j) = —0.2 (b), 
(j) = —0.5 (c). Results of the numerical simulations of Eqs. ( |8T|) with e = 0.025 and periodic 
boundary conditions. The system is 400 x 460. 
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FIG. 9: The effect of a large-scale fluctuation on a hexagonal pattern. Result of the numerical 
solution of Eq. (81) with e = 0.025, cf) = —0.2, and periodic boundary conditions. The system is 
400 X 460. 
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